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Abstract We consider polygonal Markov fields originally introduced by Arak and Sur- 
gailis [2]. Our attention is focused on fields with nodes of order two, which can he regarded 
as continuum ensembles of non-intersecting contours in the plane, sharing a number of 
features with the two-dimensional Ising model. We introduce non-homogeneous version 
of polygonal fields in anisotropic enviroment. For these fields we provide a class of new 
graphical constructions and random dynamics. These include a generalised dynamic rep- 
resentation, generalised and defective disagreement loop dynamics as well as a generalised 
contour birth and death dynamics. Next, we use these constructions as tools to obtain new 
exact results on the geometry of higher order correlations of polygonal Markov fields in 
their consistent regime. 
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1 Introduction 

The first example of a polygonal Markov field has been provided by Arak [1]. Further de- 
velopments are due to Arak & Surgailis [2, 3], Surgailis [18] and Arak, Clifford & Surgailis 
[4]. The fields with nodes of order 2, or V-shaped nodes for short, which are in the focus of 
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our attention in this paper, arise as ensembles of self- avoiding closed polygonal contours in 
the plane interacting by hard core exclusions, possibly with some further terms entering the 
Hamiltonian. Under a particular choice of the Hamiltonian the polygonal fields enjoy strik- 
ing properties including consistency (the field constructed in a subdomain D C D' C 
coincides with the restriction to D of the field constructed in D') as well as availabihty of an 
explicit formula for the partition function and some other numerical characteristics of the 
field, see [2]. Under even milder conditions one can guarantee the isometry invariance of 
the field as well as the two-dimensional germ-Markov property stating that the conditional 
behaviour of the field in an open bounded domain depends on the exterior configuration 
only through arbitrarily close neighbourhoods of the boundary, see ibidem. A particu- 
larly interesting class of processes seem to be length- and area-interacting modifications of 
consistent fields, which not unexpectedly exhibit many features analogous to those of the 
two-dimensional Ising model, including the presence of first order phase transition at low 
enough temperatures (Nicholls, [14]; Schrciber, [15]) and low-temperature phase separa- 
tion phenomenon (Schreibcr, [16]). Consistent polygonal fields and their length-interacting 
modifications have also interesting statistical applications where they are used as priors in 
Bayesian image analysis (Clifford & Nicholls, [6]; Kluszczyhski, van Lieshout & Schreiber, 
[10, 11]; van Lieshout & Schreiber, [12, 17]). 

With all the afore-mentioned similarities the polygonal fields enjoy with the Ising model 
as well as with more general lattice models used in Bayesian image analysis, one salient 
difference is that whereas the latter can be analysed using a rich mixture of combinatorial, 
analytic and geometric techniques, the overwhelming majority of crucial results for the 
polygonal Markov fields have been obtained via purely geometric methods usually going 
under the guise of graphical constructions. These were first provided by Arak & Surgailis 
[2, 3] and Arak, Clifford & Surgailis [4] in the form of the so-called dynamic representations 
for consistent polygonal fields. Later, we have introduced disagreement loop dynamics and 
contour birth and death representation for broader classes of polygonal fields, see [15]. The 
purpose of the present paper is twofold. 

• First, we introduce a class of new graphical constructions and random dynamics for 
a rather general class of non-homogeneous polygonal Markov fields with V-shaped 
nodes. These constructions include a generalised dynamic representation, generalised 
and defective disagreement loop dynamics as well as a generalised contour birth and 
death dynamics. 
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• Next, we use these constructions as tools to obtain new exact results on the geometry 
of higher order correlations of polygonal Markov fields in the consistent regime. 

We envision to use the new graphical representations in our current research in progress 
on polygonal fields, especially in context of obtaining further exact formulae for higher 
order correlations and tentatively also to obtain exact formulae for the free energy of the 
field (exact solution) at least in the rectangular case. Due to their intrinsically algorithmic 
nature the new constructions will also be used for Bayesian image segmentation purposes, 
which is our work in progress in cooperation with van Lieshout, see Schreiber & van 
Lieshout [17] for some partial developments in the particular tesselation-based set-up. 

The paper is organised as follows. In Section 2 below we introduce general non- 
homogeneous polygonal fields in non-homogeneous and anisotropic environments. This 
is done by imposing a general activity measure on the space of straight lines in the plane, 
determining how likely the presence of an edge along a given line is. The usual Ising-type 
length interaction is replaced by its anisotropic version there - the energy cost of producing 
an edge crossing hues of high activities is higher than the cost of an edge of the same length 
crossing only low activity lines. Further, in Section 3 general properties of such fields in 
their consistent regime are discussed and a non-homogeneous counterpart of the Arak & 
Surgailis [2] dynamic representation is developed. Next, in Section 4 we develop a gen- 
eralised dynamic representation for consistent polygonal fields, which allows for a fiexible 
choice of the field creation dynamics, parametrised by an increasing family of convex com- 
pact subsets of the field domain. In the subsequent Section 5 we employ the generalised 
dynamic representation to provide exact factorisation statements for higher order correla- 
tion functions of polygonal fields in consistent regime. In the next Section 6 wc introduce 
variants of the so-called disagreement loop dynamics for polygonal fields, including the 
generalised and defective disagreement loop dynamics. These are then used in Section 7 to 
establish exponential decay of higher order correlations in the particular set-up of rectan- 
gular fields. Finally, in Section 8 we introduce general contour birth and death dynamics 
for polygonal fields in low enough temperature, in the spirit of Fernandez, Ferrari & Garcia 
[7, 8, 9] and yielding a perfect simulation scheme directly from the thermodynamic limit. 



3 



2 Non-homogeneous polygonal Markov fields 



For an open bounded convex set D define the family Fd of admissible polygonal configura- 
tions on D by taking all the finite planar graphs 7 in U dD, with straight-line segments 
as edges, such that 

(PI) the edges of 7 do not intersect, 

(P2) all the interior vertices of 7 (lying in D) are of degree 2, 

(P3) all the boundary vertices of 7 (lying in dD) are of degree 1, 

(P4) no two edges of 7 are cohnear. 

In other words, 7 consists of a finite number of disjoint polygons, possibly nested and 
chopped off by the boundary. Further, for a finite collection (l) — of straight lines 

intersecting D we write Fd(Z) to denote the family of admissible configurations 7 with the 
additional properties that 7 C {j^^^ k and 7 fl is a single interval of a strictly positive 
length for each Zj, i — 1, n, possibly with some isolated points added. 

We shall also consider the subfamily Fd|0 Q Fd of empty boundary polygonal configu- 
rations in D consisting of those 7 G F^ which have no boundary vertices. For (l) as above 
we put Fb|0(O := Td{1) n Tdih,. 

For a Borel subset of A C by [[A]] we shall denote the family of all straight lines 
hitting A so that in particular [[M^]] stands for the collection of all straight hues in R^. 
We shall also write [[A]] for the family of all linear segments in hitting A. Consider a 
non- negative Borel measure A4 on the [[M^]] such that 

(Ml) M([[A]]) < 00 for all bounded Borel A C R\ 

(M2) A1([[{a;}]]) = for all x e R\ 

Below, the measure Ai will be interpreted as the activity measure on [[M^]]. Let A-^ be the 
Poisson hne process on [[M^]] with intensity measure M. and write A;^ for its restriction 
to the domain D. Further, define the Hamiltonian : Fd — > 1R+ given by 

L^(7):= Yl ^(M), TeTz,. (1) 

eeEdges(7) 
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We argue that the energy function L-^ should be regarded as an anisotropic environment- 
specific version of the length functional. Indeed, interpreting the activity A4{dl) of a line 
I hitting an edge e e Edges(7) as the hkelihood of a new edge being created along I 
intersecting and hence fracturing the edge e in 7, we note that, roughly speaking, the value 
of A^([[e]]) determines how likely the edge e is to be fractured by another edge present in 
the environment. In other words, i^'^(7) determines how difficult it is to create the whole 
graph 7 G F/) without fractures in the environment A4 - note that due to the anisotropy of 
the environment there may be graphs of a higher (lower) total edge length than 7 and yet 
of lower (higher) energy and thus easier (more difficult) to create and to keep unfractured 
due to the lack (presence) of high activity lines likely to fracture their edges. 

With the above notation, for G M further referred to as the inverse temperature 
(for mathematical convenience we admit also the unphysical negative values of inverse 
temperatures), we define the polygonal field A^'^ in D with activity measure A4 to be 
the Gibbsian modification of the process induced on F/j by A^, with the Hamiltonian L-^ 
at inverse temperature that is to say 

ryAn eG).- EE,,r„,A4.,exp(-/JL«(7)) ^' 

for all sets G C To Borel measurable with respect to, say, the usual Hausdorff distance 
topology. Note that this definition can be rewritten as 

nA-^-'' G ^7) oc exp(-/?L^(7)) J] M{dl[e]), 7 G F^, (3) 

eGEdges(7) 

where l[e] is the straight line extending e. In other words, the probability of having A"^'^ G 
d'-f is proportional to the Boltzmann factor exp{—f3L-'^{'f)) times the product of edge 
activities A4{dl[e]), e G Edges(7). Observe also that this construction should be regarded 
as a specific version of the general polygonal model given in (2.11) in [2]. The finiteness of 
the partition function 

Z^-'f'-.^E exp(-/3L-^(7)) <oo (4) 

7Grc(A^) 



for all /? G M is not difficult to verify, see (32) below. We shall also consider the empty 
boundary version of the above construction, with A'^^f arising in law as A'^''^ conditioned 
on staying within Tdm, that is to say not containing boundary vertices. It is easily seen 
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that alternatively A^^f can be defined by rewriting (2) above with Fd replaced with Fdi^ 
throughout. The corresponding partition function Z'^^^^ is always finite as well, see (31) 
below. 

In the sequel we also consider the thermodynamic limits for such fields, which take 
their values in the space F := r^2 of whole-plane admissible configurations, with obvious 
meaning of this notation. Note in this context that r^2 = r]R2|0. 

As a particular case of this general construction we shall consider rectangular fields. 
Slightly abusing the terminology, by a rectangular field we shall mean a polygonal field 
given by (2) and (3) whose activity measure Ai concentrates on translates of two non- 
parallel straight lines - in other words, its edges can follow precisely two different directions. 
Clearly, upon a non-singular affinc transform such a rectangular field can be made into a 
process where all angles between edges are straight, whence the name. Assuming with no 
loss of generality that the two admissible directions for the edges of the field are parallel 
to the coordinate axes, we can further transform each rectangular field into the standard 
rectangular field where the measure A4 is given by Ai{dl) being dh if / is a parallel translate 
of a coordinate axe by distance h, /i e M, and otherwise. Indeed, putting 



F / M((0,x]x{0}), iix>0, 

J'^K-^) ■ \ -M{{x,0]x {0}), otherwise 



and 



A^({0}x (0,1/]), ify>0, 
-M{{0} X {y,0]), otherwise, 

it is readily verified that the map 3 {x, y) ^ {F^{x), Fi{y)) does define such a transfor- 
mation. Note that both and are finite by (Ml) and continuous by (M2). Although 
in general this transformation is not a bijection, it is easily seen to be bijective on the sup- 
port of A^, with the inverse defined by {s,t) i— > (sup{x, F^(x) < s},sup{y, -Fj(y) < 0)- 
Even though in general we admit all /3 e IR, in this paper we shall mainly study two 
principal temperature regions. The first one is the consistent regime corresponding to P — 1 
and introduced in Section 3 - we shall argue that this particular choice of temperature 
parameter places us in the context of a non-homogeneous version of Arak-Surgaihs [2] 
construction for the so-called consistent polygonal fields, thus ensuring the availability of 
an appropriate dynamic representation for our process in terms of equilibrium evolution 
of one-dimensional particle systems tracing the boundaries of the field in two-dimensional 
space-time, as discussed in detail in Subsection 3.1 below, see also [2]. The afore-mentioned 
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consistency property, arising as a consequence of the dynamic representation and further 
discussed in Subsection 3.2, means here that for each open bounded convex D C D' c 
the field A'^'^ :— A^'^ coincides in law with A-^, :— A^,'^ fl D, thus allowing for a 
direct construction of the infinite volume process (thermodynamic limit) A^^ :— A^^''^ 
on the whole R^. A number of further properties can be concluded from the dynamic 
representation for (3 — 1, including the explicit knowledge of the partition function Z^'^ 
as made precise in Theorem 2 below. The consistent fields are the main object of study in 
Sections 3, 4, 5 and 7, and play a crucial role also in Section 6. 

The second regime in the focus of our interest is the low temperature region (large 
positive f3) where long range point-to-point correlations arc present, giving rise to the 
spontaneous magnetisation phenomenon, see Corollary 6 in Section 8 and the discussion 
below it. All disagreement loop dynamics developed in Section 6 of this paper admit their 
versions working for the low temperature region. Moreover, in Section 8 wc provide another 
important tool characterising the behaviour of our model in this regime, namely a contour 
birth and death graphical construction in the spirit of Fernandez, Ferrari and Garcia [7, 8, 9] 
put together with our random walk representation for the so-called free contour measure 
specific to the model. Applying this graphical construction we are able to establish the 
existence of the thermodynamic limit for A"^'^ as well as certain mixing-type results for 
general activity measures satisfying (27) in analogy with Section 4 in [15], see Corollary 5. 

As a by-product of our considerations both in the consistent and low temperature 
regimes we obtain efficient exact simulation algorithms for A^^''^, both in finite windows 
and directly from the thermodynamic limit. 

3 The consistent regime 

As mentioned above, we show here that the choice (3 — 1 puts us in the so-called consistent 
regime of the polygonal field, where it enjoys remarkable features arising as consequences of 
a (non-homogeneous) version of the dynamic representation developed by Arak & Surgailis 
[2], see Section 4 there. In Subsection 3.1 below, we give the non-homogeneous dynamic 
representation and further in Subsection 3.2 we discuss its consequences. 
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3.1 Consistent fields and their dynamic representation 

To describe the dynamic representation, we interpret the open convex domain D as a set of 
space-time points (t, y) G D, with t referred to as the time coordinate and with y standing 
for the spatial coordinate of a particle at the time t. In this language, a straight line segment 
in D stands for a piece of the space-time trajectory of a freely moving particle. In the sequel 
by a particle we shall understand a function t ^ yt assigning to each time moment t the 
spatial location yt of the particle at the time t. Moreover, by the trajectory of a particle we 
shall mean the graph of this function, i.e. the set of points (t, yt). The particle trajectories 
considered in this paper are piecewise linear and hence almost everywhere differentiable. 
Clearly, in the language being introduced, the derivative of yt at t corresponds to the current 
velocity of the particle at the time t. Below, for notational convenience we shall omit the 
subscript t writing y instead of yt- For a straight line I non-parallel to the spatial axis and 
crossing the domain D we define in the obvious way its entry point to in(Z, D) e dD 
and its exit point out(Z, D) e dD. 

On we construct the hirth measure {{M)) given for each Borel ^4 C by 

{{M)){A) := ]-M X M{{{li, h): hnkC A}) = Ecard{{Zi, k} C A^, h + hnhc A}. 

^ (5) 

Note in this context that M x M{{{1,1), I E [[K^]]}) = by (M2). Likewise, on dD we 
construct the boundary birth measure {{M.; dD)) given for each Borel B C dD by 

{{M;dD)){B) := M{{1, in(/, D) e 5}) = Ecard{/ G A"^, m{l,D)eB}. (6) 

We choose the space-time birth coordinates for the new particles according to a Poisson 
point process in D with intensity measure ((A^)) (interior birth sites) superposed with a 
Poisson point process on the boundary (boundary birth sites) with the intensity measure 
{{M; dD)). Each interior birth site x E D emits two particles, moving with initial velocities 
v' and v" chosen according to the joint distribution 

nMrj f J //N ^ X -^({(^i^ k), knke dx,v[k] e dv', vlk] e dv"}) 
' ) ^= 2{{M)){dx) 

where v[l] is the velocity of a particle whose space-time trajectory is represented by the 
straight line /. Note that this is equivalent to choosing the directions of the straight lines 
representing the space-time trajectories of the emitted particles according to the distribu- 
tion of the typical angle between two lines of A-^ at x, that is to say the angle arising by 
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conditioning A-^ on containing two lines intersecting at dx, see also Sections 3 and 4 in [2] 
and the references therein. In other words, (7) implies that the interior birth at li fl Z2 in D 
of two particles moving along straight lines h and I2 respectively happens with probability 
A4{dli)A4{dl2). Each boundary birth site x e dD yields one particle with initial speed v 
determined according to the distribution 9^''^'{dv) identified by requiring that the direc- 
tion of the line entering D at x and representing the space-time trajectory of the emitted 
particle be chosen according to the distribution of a straight hne I e A-^ conditioned on 
the event {x = in(/, D)}, that is to say 

nM■,^f,^ . m{l,D) G dx,v[l] e dv}) 
^^"^-^ {{M;dD))idx) ■ 

Again, this means that the boundary birth at I fl dD of a particle moving along a straight 
hne I occurs with probability M.{dl). 

All the particles evolve independently in time according to the following rules. 

(El) Between the critical moments listed below each particle moves freely with constant 
velocity so that dy — vdt, 

(E2) When a particle touches the boundary dD, it dies, 

(E3) In case of a collision of two particles (equal spatial coordinates y at some moment 
t with {t, y) E D), both of them die, 

(E4) The time evolution of the velocity Vf of an individual particle is given by the following 
pure-jump inhomogeneous Markov process: denoting by It the straight line extending 
the present segment of the space-time trajectory of the particle at time t so that 
Vt — v[lt], we have 

F{lt+dt e dl' \lt = l) = l{i,^^^iit,t+,t^}{{M)){dl') 
with standing for the segment along I between time moments t and t -\- dt. 

The evolution rule (E4) can be interpreted as follows: in the course of its movement along 
a short line segment a particle turns and starts moving along another straight line 

I' e ]] with probabihty {{M.)){dl') which agrees with the interpretation of as 

the line activity measure as discussed above. 
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The following theorem shows that the polygonal field obtained in the course of the 
above dynamic construction is in fact A'^ — A^'^, in analogy with the results of Arak & 
Surgailis [2], Section 4, for homogeneous fields. 

Theorem 1 The polygonal field traced by the particle system constructed above coincides 
in distribution with A"^ — A^'^ . Moreover, we have 

Z^-'' = eM{{M)){D)). (9) 

Proof We pick some 7 G Fd and calculate the probability that the outcome of the above 
dynamic construction falls into dj. To this end, we note that: 

• Each edge e e Edges(7) with initial vertex (lower time coordinate) lying on dD 
contributes to the considered probabihty the factor A4{dl[e\) (boundary birth of a 
particle tracing the edge) times exp(— Al([[e]])) (no velocity updates along e). 

• Each of the two edges 61,62 G Edges(7) steming from a common interior birth site 
l[ei] n ^[62] yields the factor Ai{dl[ei]), i = 1,2, (coming from the birth probability) 
times exp(— A^([[ei]])) (no velocity updates along Cj), 

• Each of the edges e G Edges(7) arising due to a velocity update of a particle yields 
the factor A4{dl[e]) (velocity update probability) times exp{—A4{dl[e])) (no velocity 
updates along e), 

• The absence of interior birth sites in D \ 7 yields the factor exp(— ((A1))(D)), 

• Finally, the absence of boundary birth sites at dD \ 7 yields the additional factor 
exp{-{{M;dD)){D)) = exp{-MiD)). 

Putting these observations together we conclude that the probabihty element of 7 being 
traced by the particle system is 

exp(-L^^(7)) ne.Ed,es(.)-M(^/N) ^ exp(-L^(7)) ^ 
exp{{{M)){D)) exp{M{D)) exp(((A<))(D)) ^ 

where the event {7 G roi^)'^} is easily seen to correspond to the situation where the 
collection of edge-extending lines {l[e], e G Edgcs(7)} coincides with the collection of lines 
determined by A-^. This shows that the field traced by the particles coincides in law with 
A'^ and that (9) holds, thus completing the proof of the theorem. □ 
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3.2 Conclusions from the dynamic representation 

The crucial property of polygonal fields for (3 = lis their consistency stated in Theorem 
2 below and established using the dynamic representation from Theorem 1 in analogy with 
Arak & Surgailis [2]. This is where the name consistent regime comes from. Another 
important feature of is the explicit knowledge of its linear sections, see Theorem 2, 
which can be further used to obtain information about the first and second order correlation 
structure of the field, see Theorem 4 and Corollaries 1, 2 and 4. 

Theorem 2 The polygonal field enjoys the following properties 

Consistency: For bounded open convex £)' C D C the field A'^ fl D' coincides in law 
with A^i. This allows us to construct the whole plane eoctension of the process A^ 
such that for each bounded open convex D C the field A'^ coincides in law with 

A-^r\D. 

Linecir sections: For a straight line I in the collection of intersection points and 
intersection directions of I with the edges of the polygonal field A''^ coincide in law 
with the corresponding collection for the Poisson line process A^. 

Proof To establish the Consistency property, choose a bounded open convex set D C 
R^ and a straight fine I intersecting D and define D' to be the set of points of D lying to 
the left from I (lower time coordinates). Clearly then, from the dynamic representation 
Theorem 1 we conclude the Consistency statement for the so chosen D and D' . Noting 
that the dynamic representation is equally available upon rotating the space-time coordi- 
nate system we see that the Consistency holds as well upon cutting off the part of the 
set D lying to the left of I. This means however that the consistency holds upon cutting off 
pieces of the original set with arbitrary straight lines - a repetitive use of this procedure 
and a possible passage to the limit allows us to carve from D its arbitrary convex subset. 
This proves the Consistency claim. To estabfish the Lineeir sections statement pick 
a straight line I and choose the space-time coordinate system so that / coincides with its 
spatial axis. The Linear sections claim follows now from the form of the boundary birth 
mechanism in the dynamic representation in view of the Consistency property. □ 
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4 Generalised dynamic representation for consistent 
fields 



The dynamic construction of consistent polygonal fields, originating in the homogeneous 
set-up from Arak & Surgailis [2] and discussed in Subsection 3 above, can be regarded as 
revealing increasing portions of the polygonal field in the course of the time flow. Under 
this interpretation, the portion of a polygonal field in a bounded open convex domain 
D uncovered by time t is precisely its intersection with Dt — D n {—oo,t\ x R+. The 
idea underlying our generalised version of dynamic representation developed in the present 
section is to replace the above family Dt by some other time-increasing family of subsets 
of D, also denoted in the sequel, eventually covering the whole D, and to try to provide 
a natural construction of the polygonal field being gradually uncovered on the growing 
domain Df in the course of the time flow. We shall always assume that Df be convex for 
otherwise we would have to deal with situations where two or more disconnected parts of an 
edge of the fleld have been revealed which leads to unwanted and cumbersome dependencies 
along the segments connecting these parts. Taking this into account and having formal 
convenience in mind we impose the following natural assumptions on D and Dt, t & [0,1], 

(Dl) (Dt)te[o,i] is a strictly increasing family of compact convex subsets oi D — DL) dD. 

(D2) Dq is a single point x in D — DU dD. 

(D3) Dl coincides with D. 

(D4) The family {Dt)te[o,i] enjoys the property that M{{1, ^telo,i] card(Zn9A) > 2}) = 0, 
in particular D itself satisfies A4{{1, card(Z fl dD) > 2}) = 0. 

(D5) Dt is continuous in the usual Hausdorff metric on compacts. 

Under these conditions, for Al-almost each I e [[£)]] the intersection I fl D^-, consists of 
precisely one point A(/), where ti — mi{t e [0, 1], Dt D I ^ 0}. The point A(Z) will be 
referred to as the anchor point for I, this induces the anchor mapping A : [[D]] D. 
Consider now the following dynamics in time t G [0, 1], with all updates given by the rules 
below performed independently of each other. 

(GErlnitialise) Begin with empty field at time 0, 
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(GEl) Between critical moments listed below, during the time interval [t,t + dt] the 
unfolding field edges in Dt reaching dDf extend straight to Dt+dt \ Dt, 

(GE2) When a field edge hits the boundary dD, it stops growing in this direction (recall 
that A^-almost everywhere the intersection of a line with dD consists of at most two 
points) , 

(GE3) When two unfolding field edges intersect in Dt+^t \ Dt, they are not extended any 
further beyond the intersection point (stop growing in the direction marked by the 
interesction point), 

(GE4) A field edge extending along I e [[Dt]\ updates its direction during [t,t + dt] and 
starts unfolding along I' e extending away from the anchor point A(Z'), 

with probability M{dl'), where ;[*'*+<^*l :— I n {Dt+dt \ Dt). Directional updates of this 
type are all performed independently, 

(GEiLineBirth) Whenever the anchor point A{1) of a line I falls into Dt+dt \ Dt, the 
line I is born at the time t at its anchor point with probability M.{dl), whereupon it 
begins extending in both directions with the growth of Dt (recall that I is Al-almost 
always tangential to dDt here) , 

(GE:VertexBirth) For each intersection point of lines li and I2 falling into Dt+dt \ Dt, 
the pair of field lines li and I2 is born at li fl I2 with probability Ai{dli)A4{dl2), 
whereupon both lines begin unfolding in the directions away from their respective 
anchor points A(/i) and A{l2). 

Observe that the evolution rule (GErVertexBirth) means simply that pairs of fines are 
born at birth sites distributed according to a Poisson point process in D with intensity 
measure {{M)), in analogy to the standard dynamic representation in Subsection 3.1. The 

essential difference though is that here the pairs of lines emitted from a birth site extend 
away from their respective anchor points rather than always in the direction determined 
by a single time axis. Further, the line birth events given by (GE:LineBirth) replace 
the boundary birth events in the standard dynamic representation. It is also worth noting 
that if we choose the family Dt so that Dt := D (1 (—00, (1 — t)x„ii^ + temax] x where 
Xmin and Xmax are the minimal and maximal x-coordinates of a point in D (assume that 
D contains exactly one point with x-coordinate Xmin and that M. assigns zero mass to 
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the set of vertical lines), the generalised dynamic representation (GE) coincides with the 
standard Arak & Surgailis one determined by rules (El-4) in Subsection 3.1 above, and 
we have A{1) — m{l, D). 

In analogy with the corresponding result for the usual dynamic construction, as estab- 
lished in the proof of Theorem 1, we show that the field resulting from the above (GE) 
construction does coincide in law with A'^. 

Theorem 3 The random contour ensemble resulting from the above construction (GE) 
coincides in law with A^. 

Proof The proof is very similar to that of Theorem 1. We pick some 7 G and calculate 
the probability that the outcome of the above dynamic construction falls into ^7. To this 
end, we note that: 

• Each edge e G Edgcs(7) containing the anchor point A(/[e]) and hence resulting from 
a line birth event due to the rule (GE:LineBirth), contributes to the considered 
probability the factor M.{dl[e]) (line birth probability for l[e]) times exp(— A1([[e]])) 
(no directional updates along e), 

• Each of the two edges ei, 62 G Edges (7) steming from a common interior birth vertex 
Z[ei] n ^[62] yields the factor M.{dl[ei\), i — 1, 2, (coming from the vertex birth prob- 
ability due to the rule (GEiVertexBirth)) times exp(— A^([[ei]])) (no directional 
updates along e^), 

• Each of the edges e G Edges(7) arising due to a directional update in (GE4) yields 
the factor A4{dl[e]) (directional update probability) times exp{—J\4{dl[e])) (no direc- 
tional updates along e), 

• The absence of interior birth sites in D \ 7 yields the factor of exp(— ((A1))(D)), 

• Finally, the absence of line birth events for all lines in [[D]\ except for the finite 
collection {Z[e], e G Edges(7)} yields the additional factor exp(— Al(-D)). 

Putting these observations together we conclude that the probability element of 7 resulting 
from the generalised construction above is 

exp(-L^(7)) UeeEd,esMM{dl[e]) _ exp{-L^{j)) 



eM{{M)){D)) eMM{D)) exp(((A^))(D)) 



rfP(7 G r^(A^)) 
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and thus the field obtained by this construction coincides in law with as required. 
This completes the proof of the theorem. □ 

5 Correlations in the consistent regime 

In the present section we use the dynamic construction and its generalised version to 
describe the correlation structure of the consistent field A^. Due to the polygonal nature 
of the considered field the natural object of our interest are the edge correlations 

a^[dh, xi,... ; dlk, Xk] := P {yt=i^eemges{A^) ^hixi) e e, /[e] e dk) , (10) 

where h, . . . ,lk are straight lines and tt;. is the orthogonal projection on 1^. In almost all 
cases below we shall be interested in correlations with Xi & k, in which case a-'^[dli, xi, . . . ; 
dlk,Xk] can be interpreted as the probability element that the polygonal field A'^ passes 
through points Xi in the directions determined by the respective lines Zj, i = 1, . . . ,k. 
For general Xi, not necessarily lying on li, the /c-fold correlation (T-^[(i/i, Xi; . . . -jdlkyXk] 
is the probability that the polygonal field passes through points 7ri.(a;j) in the directions 
determined by the respective lines /j, i = 1, . . . ,k. 

5.1 First and second order edge correlations 

The first and second order edge correlations are easily determined using the Lineair sec- 
tions statement of Theorem 2 concluded from the standard dynamic representation. 

Theorem 4 We have for x e I 

a^[dl,x]^M{dl) 

and, for Xi E li, X2 E I2, 

a^[dli,xi;dl2,X2] = 
where X1X2 is the segment joining Xi to X2- 



M{dh)M{dl2), tfh^k, 
eM-'^M{\[x^]]))M{dl), ifh = I2, 



Proof The statement for a-^[dl, x] and a-'^[dli,xi] dl2, X2], /i 7^ /2, is a direct consequence 
of the Linear sections part of Theorem 2 and of the properties of the Poisson line process 
A-^. To find [dl, Xi;dl, X2] for xi,X2 E I note that in order to have an edge e of the field 
passing through both xi and X2 along I, we have to ensure that 
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• There is e e Edges(^-^) such that l[e\ e dl and xi, which happens with probabihty 
M.{dl) in view of Linear sections in Theorem 2, 

• There are no other edges of the field crossing e between Xi and X2, which happens 
with probabihty exp(—2A^ ([[0:1X2]])). Indeed, to see it choose the spatial axis in the 
standard dynamic representation in Subsection 3.1 very close to I (not exactly I to 
keep the velocity of the particle tracing e finite though very large) and observe that 
e can be crossed by 

— space-time trajectories of particles coming from the past, which happens with 
probability 1 — exp(— A1([[aJilr2]]))[l + o(l)] by the Linear sections property, 

— edges arising due to velocity updates along X1X2, which happens with probability 
1 - exp(-A^([[xIx^]]))[l + 0(1)] by the dynamic rule (E4). 

Since the velocity update events in (E4) arc independent of the past particle config- 
uration, letting the spatial axis approach / we obtain the required conclusion. 

Combining the factors listed above we obtain the required formula for a-'^[dl,xi;dl,X2\, 
thus completing the proof of the Theorem. □ 



5.2 Higher order edge correlations for general activity measures 

To describe the higher order correlation structure of the field we need the full power of 
the more flexible generalised dynamic representation. Consider a collection {li,xi), {I2, X2), 
. . . , {lk,Xk), k > 1, oi pairwise different lines k and points Xi with Xj G k and Xi ^ Ij for 
j 7^ i. Such collections are said to be in general position below and such assumption will 
be imposed an all collections {li,Xi) considered in this subsection, often without a further 
mention. Also, throughout this subsection we always assume for formal convenience that 
all contained in a bounded open convex set playing the usual role 

of the field domain. Clearly, the particular choice of D is irrelevant due to the consistency 
of the field 

We say that the edge correlations of the field A'^ factorise on a collection (Zj, Xj)^^^ if 
a^[dli, xi, . . . ;dlk, Xk] coincides with the product Y[i=i o''^[dli, Xi\. Moreover, we say that 
the collective factorisation of correlations holds for a family {li,Xi)^^i if the correlations 
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factorise for (Zj, Xi)'^^^ and all its sub-collections. In view of Theorem 4 (Zj, Xi)^^i collectively 
factorises iff 

m 

a^[dk„Xi,;...;dk^,XiJ^l[M(dk.) (11) 
for all sub-collections (^i^ , a^i^)^!- 

General collective factorisation problem and precedence graphs Our main ob- 
jective in the present subsection is to find general conditions characterising collections 
{l,x) in general position enjoying the collective factorisation property. To see what may 
be plausible answers to this question let us make the following basic observations very 
helpful in interpreting the geometry of the higher order edge correlations. For a collection 
(l, x) = {li,Xi; . . . ;lk, Xk), Xi e in general position by T(I, x) we shall mean the family 
of admissible polygonal configurations 7 in the plane consisting of k edges ei, . . . , such 
that ei lies on k and contains for alH = 1, . . . ,k. Looking at the structure of the family 
r{l,x) is closely related to the following natural problem which marks its presence in vari- 
ous domains ranging from kid games to studies on human and computer vision, see [5] and 
the references therein: given the collection {l,x) = {li,Xi)f^^ draw a family of closed curves 
(here a polygonal configuration) such that each point Xi lies on one of the curves and the 
direction of the curve at Xi is determined by /j. Clearly, the solution to this problem is 
non-unique in general, yet if we require in addition that the configuration we draw belong 
to r{l,x), that is to say we may only draw over the lines li and on each of these lines we 
have to draw precisely one segment of non-zero length, then it is often the case that r(/,a;) 
is a singleton. Now, assume that all points Xi as well as the intersection points yij for k 
and Ij lie very close to each other, say they are all contained in a disk of very small radius 
r, then it follows from (2) and (3) combined with the consistency property of the field that 

a^[dh,xi;... ]dlk,Xk] = N{l,x)M{dh) . . . M{dlk){l + Or{l)), (12) 

where N(l, x) is the cardinahty of V{1, x). Indeed, to see it note that the Boltzmann weight 
exp(— L-^(-)) tends to 1 as r — > 0. Thus, the factorisation holds in small r asymptotics 
precisely when T{1, x) is a singleton. As will be discussed in the sequel, there are many 
collections for which N(l,x) = or N{l,x) > 1. In particular, it is not the case that 
the factorisation holds for all collections in general position. Neither can it be hoped 
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though that the formula (12) holds in general non- asymptotic regime: as we shall show in 
Theorem 11 in Section 7 in the particular case of rectangular fields, when the distances 
between Xj's get large, the correlations converge exponentially fast to the product, that 
is to say the polygonal field exhibits exponential decay of dependencies. To conclude 
these considerations, we say that {li,Xi)'^^i enjoys collective factorisation on all scales iff 
(a/j, axi)^^^ factorises collectively for each a > 0, where by {al, ax) we understand the 
re-scaled version of {l,x) with scaling factor a. Then the above discussion shows that 

Lemma 1 For a collection {l,x) — {li,Xi)^^i in general position a necessary condition to 
enjoy collective factorisation on all scales is that N{l,x) — 1 and N{l',x') — 1 for all 
non-empty sub-collections {I', x') of {I, x). 

To reformulate this condition in more tangible terms, we build for each family (/j, Xi)^^i its 
precedence graph Xi; . . . ; Z^, Xk] as follows. We split each at Xi into two half- lines both 
oriented in the directions away from Xi. This creates a directed graph Q[li,Xi; . . . ; lk,Xk] 
whose vertices are the generating points Xi and the intersection points yij of respective pairs 
of lines /«, Ij. As we shall see below, the name precedence graph comes from its relationship 
to the order in which the points of Ijf^]^ k are revealed in the course of a suitable instance 
of the generalised graphical construction. For now, we claim that 

Lemma 2 For a collection {l,x) = in general position the following conditions 

are equivalent 

1. N{l,x) — 1 and N{l',x') — 1 for all non-empty sub-collections of {l,x), 

2. The precedence graph Q[li,Xi] . . . ; x^] is acyclic. 

Proof To show that the first condition implies the second one note that if the precedence 
graph contains a cycle of length m built by (/j^, Xj^; . . . ; li^,x.i^), then for the sub- collection 
{l',x') = (/ij , ) we have N{l',x') = if m is odd and N{l',x') = 2 if m is even. 

To prove the inverse implication assume that the precedence graph Q[li,Xi; . . . ; lk,Xk] 
contains no cycles and observe that this induces a partial ordering on IJi=i ^* which the 
Xj's are minimal points, each Xi first on its corresponding line li and with the remaining 
points on /j ordered according to the natural linear orderings directed away from Xi on the 
two half- lines. This will be referred to as the structural ordering for the collection (/, , x'^jLi 
in the sequel. Consider the following incremental construction of a graph belonging to 
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r{l, x), where at each step we obtain a collection of k segments Li C. k, i — 1, . . . ,k with 
non-intersecting interiors but possibly sharing endpoints and possibly degenerated to Xj's, 
eventually yielding the entire graph under construction. 

1. begin with ij :— {xi}, i — 1, . . . ,k, 

2. choose an intersection point Uij which: 

(a) belongs to neither tj nor ij and hence to neither of the remaining segments li. 

(b) enjoys the property that extending both ij and Lj to contain y^j creates no T- 
shaped or X-shaped nodes (vertices of order three or four respectively), which 
is equivalent to both ^ and Lj having their endpoints pointing towards yij not 
shared with any other segment i/. (Such endpoints will be referred to as loose 
ends below. Note that by definition of our construction procedure the only 
possible loose end on l^ is Xi.) 

(c) is minimal with the above two properties, 
with possible ties broken in an arbitrary way, 

3. extend both segments Li and Lj to contain y^j - this action will be referred to as 
adding i/ij to the graph for short in the sequel of this argument, 

4. return to 2. unless no more UijS can be added, 

5. whenever a segment li of the constructed graph ends with a loose node, that is to 
say a vertex of order one, extend this segment to the half-line in the direction of the 
node. 

In other words, we initialise our graph under construction with the set of generating points 
{xi, . . . ,Xk} whereupon we let it grow along the lines li, adding subsequent intersection 
points Uij in the order determined by the structural ordering whenever this does not 
violate the usual constraints imposed on a polygonal configuration, and discarding those 
yjj's whose addition would violate these constraints. It is easily seen that this procedure 
yields in a finite number of steps a graph belonging to r{l, x), thus in particular r{l, x) ^ 
and N(1,x) > 1. To show that N(l,x) < 1 note that in fact all graphs in r(l,x) can be 
obtained by the above procedure. Indeed, this is seen inductively. First, all x^'s have to 
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belong to such a graph. Second, whenever a graph in r{l,x) contains some configuration 
of tj's as its (sub) segments and some of these tj's have loose ends Xj's, each of such tj's has 
to be extended, either up to intersection with the extension of another segment or to the 
entire half-line if such an intersection is not available. This is because no loose ends can 
be present in the final graph and because we cannot hope that the potential loose end Xi 
on Li could be possibly reached by another edge ij, j i, since Xi ^ Ij, j ^ i. This means 
that at least one intersection point yij as in 2.(a,b,c) has to be added to the graph as in 2. 
Consequently, two different graphs in r(/,a;) could only arise in the course of two distinct 
instances of the above procedure, with different choices of minimal points yij in step 2. 
We argue that the so constructed graphs necessarily coincide though. To see it let us make 
the following observations 

• If at a certain stage of the construction some intersection point i/ij becomes admissible 
for 2.(a,b,c), there is no way to make it inadmissible whatever choices be made in 
the sequel of the construction save for adding i/ij to the graph. Indeed, to make yij 
inadmissible in the sequel of the construction we would have to close either tj or ij 
by intersecting its extension with an extension of some other li before tj or ij hits 
yij, which is not possible by the minimality of ytj in 2.(a,b,c). 

• If at a certain stage of the construction some intersection point yij becomes admissible 
for 2.(a,b,c), there is no way to make it inadmissible by adding some other intersection 
points in the prequel of the construction in agreement with the rule 2. (without 
removing those already present though). Indeed, to make yij inadmissible in this 
way we would have to have added some other intersection point y on tj or lj prior to 
the considered construction stage, but if this were possible this point y would have 
to reach its admissibility before yij and in view the previous observation there would 
be no way to make it inadmissible before yij reaches its admissibility, which would 
contradict the minimality of y^j in 2.(a,b,c) at the moment of becoming admissible. 

Consequently, regardless of any particular sequence of choices in a given instance of our 
construction, the points admissible for 2.(a,b,c) in the very first step of the construction 
(call them the first generation) eventually have to be added in view of the first observation 
above. Next, by the first and second observation above, all points which become admissible 
upon adding just the first generation (call them the second generation) also have to be 
added at some stage of the construction regardless of the particular sequence of choices 
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made. Proceeding further this way we conclude inductively that the so-defined generations 
of all orders will eventually be added to the graph. However such generation-wise addition 
of intersection points (plus possible extensions to half-lines at the end of the procedure) is 
also an instance of our incremental construction and it yields a valid graph belonging to 
T(1,x) which, in view of the above discussion, enjoys the property of being contained in 
any other graph obtained in any instance of our construction. It remains to observe that a 
valid graph in T(J, x) cannot be further extended by our construction. All this means that 
the results of all possible instances of our construction coincide and the order in which the 
points are considered in 2. is irrelevant. This way, we have shown that N{l,x) = 1 as 
required. □ 

We ask whether the necessary precedence graph acyclicity condition stated in combined 
Lemmas 1 and 2 is also a sufficient condition for collective factorisation of correlations (on 
all scales). Below we are going to show that the answer to this question is positive for 
rectangular fields, see Subsection 5.3 and Theorem 6 there. For general fields we were 
only able to establish collective edge factorisation under a somewhat stronger sufficient 
condition though, see Theorem 5, and we do not know at present if this condition can be 
weakened. 

General sufRcient condition for collective factorisation We proceed with the gen- 
eral case first. To this end, we note that the generalised dynamic construction with its 
(GE:LineBirth) rule allows us to conclude: 

Lemma 3 Assume that {li,Xi)f^i is a collection of pairwise distinct lines and points in 
general position and such that h does not hit the convex hull conv({a;2, ... ,Xk})- Then 

a-^[dli,xi] . . . •,dlk,Xk\ = M{dli)a-^[dl2,X2; . . . •,dlk,Xk\- 

Indeed, under the assumptions of the corollary, by standard geometry it is always possible 
to construct an increasing family (Dt), t G [0, 1], of convex compacts satisfying rules (Dl- 
5) as in Section 4 and such that Dt covers the whole of conv({a;2, . . . , x^}) before hitting li 
and that A(/i) = xi, that is to say xi is the first point of li hit by Dt. This yields the required 
statement by the consistency of the field and by the (GE:LineBirth) rule. This way of 
thinking suggests a natural sufficient condition for collective factorisation to hold under 
general activity measures, namely that there exist an increasing family (Df), t e [0, 1], of 
convex compacts as in Section 4 and such that Xi — A{li) for all 1 < i < k, that is to say Xi 
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is the point at which li is first hit by Df. However, the so formulated collective factorisation 
condition is rather untractable, therefore we look for its equivalent reformulation in more 
tangible terms. To this end we augmentthe precedence graph xi, . . . ; Ik, Xk] as follows. 
For each two lines li, Ij intersecting at some yij we note that Zj and Ij divide the plane into 
four regions, one with both Xi and Xj on its boundary, two with either Xi or Xj, and finally 
one, termed the trap region, with neither. Now, if some Xm falls into such a trap region, 
it is easily seen that for each family Dt such that A(/j) = Xj, i = 1, . . . ,k, we must have 
yij hit by Df prior to x^- Indeed, if at a certain time Dt contains x^, then if by that time 
it contains either of the points Xj or Xj, it has to intersect both /j and Ij and thus contain 
both Xi, Xj and hence also yij by convexity. On the other hand if Dt hits Xm before hitting 
any of Xi and Xj then we cannot have simultaneously A(/j) = Xi and A(/j) = Xj. Taking 
this into account we add in Q[li,xi; . . . ;lk, Xk] a directed trap edge from y^j to Xm for each 
i,j,m as above. Denote the resulting augmented directed graph by . . . -jl^jXt]. 

It is important to observe that, as follows by its construction and the above discussion, 
the orientation of edges in this graph indicates the order in which its vertices are hit by 
the sought for increasing family D^, should it exist. Consequently, the acyclicity of the 
augmented precedence graph is a necessary condition for the existence of such Df. In the 
proof of the following crucial theorem we show that it is also a sufficient condition. 

Theorem 5 Assume that {li,Xi)i^i is a collection of pairwise distinct lines and points in 
general position such that the augmented precedence graph Q^[li,Xi; . . . ; lk,Xk] is acyclic. 
Then {li,Xi)^^^ admits collective factorisation of correlations on all scales. 

Proof In context of the discussion above and in view of the generalised graphical con- 
struction and its (GE:LineBirth) rule as used in Lemma 3, it suffices to show that there 
exists an increasing family Dt such that A(/j) = Xi for alH = 1, . . . ,k. Now, to establish 
the existence of such a family it is enough to know that having Q^[li, Xi, . . . ^Ik, Xk] acyclic 
implies the existence of a permutation (si)*L]^ of indices such that /^..^-^ does not hit the 
convex hull conv({xs^, ... ,Xs^}), i = 1, . . . ,k — l. Indeed, the family Dt is then easily con- 
structed by induction in /c : it starts growing from Xg^ whereupon it hits the consecutive 
points Xs2,Xs3, .... and since Ig^^-^ is disjoint with conv({a;s^, . . . ,Xsi^}), the family Dt can 
be chosen so that it does not hit /^^^^ until it reaches all previous lines /g., i < k, and 
so that A(/sJ = Xg^ for i < /c (by induction hypothesis) and then A(/5j.^J = Xg,^^-^ (again 
by disjointness of lsi^_^_^ with the convex hull of preceding Xg.'s). This can be equivalently 
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interpreted as iterative application of Lemma 3. 

To proceed, we use induction in k to show that whenever ^+[. . . ] is acychc, the required 
permutation ensuring the disjointness of hnes with convex hulls of sets of preceding points 
exists. To this end, observe first that our statement triviahses ioi k — 1. Next, choose a 
collection {li, Xi; . . . ; Ik, Xk) and note that if there exists a line in the collection which does 
not hit the convex hull generated by all the remaining points then by Lemma 3 we can 
resort to the inductive hypothesis for the collection with the hne removed. Consequently, 
we only have to show that if the collection . . . ; Iki^k) is such that each line li hits 

the convex hull conv({xj, j ^ i}) then Q^[li,xi] . . . ^Ik^Xk] contains a cycle. Again re- 
sorting to inductive argument if needed, we may assume without loss of generality that 
(Zi, Xi; . . . ; Xfc) is minimal with this property, that is to say it does not contain a proper 
subcollection such that each its line hits the convex hull generated by the remaining points. 
This minimality assumption implies that all Xj's are extreme points (vertices) of the con- 
vex hull conv({xi, ... , Xfc}), for otherwise removing a non-extreme point we would obtain a 
subcollection enjoying the considered property. For formal convenience we let the vertices 
Xi be ordered clockwise along the boundary of the convex hull and we interpret the indices 
modulo /c, that is to say Xk+i — xi etc. By minimality, should we remove a vertex Xi from 
the collection, there exists Xj, j ^ i, such that Ij does not hit conv({x^, m ^ How- 
ever, by the assumed properties, the only possible choices for such j are j — i+e for e — +1 
or £ = — 1 and the line Ij has to cross the segment XiXi-e- Consequently, each vertex Xi is 
cut offirom the convex hull conv({xi, . . . , Xk}) by a line passing through its neighbouring 
vertex and crossing the opposite neighbouring edge. The choices of e for different vertices 
are not independent, because if we choose to use to cut off Xj+i then has to be cut off 
by Zj_2. Thus, if the number k of vertices is odd, this is easily checked to imply that either 
for alH = 1, . . . ,k the line li passes through the segment XjXj+i or for alH = 1, . . . ,k the 
line li passes through the segment XiXiZi. In both cases this generates a cycle of order k in 
Q'^[li,Xi; . . . ; Ik, Xk\- On the other hand, if k is even, apart from the above two options gen- 
erating a cycle of order k we have another possibility where the collection of vertices splits 
into k/2 pairs of neighbours Xi, Xj+i such that k crosses Xi^^iXj^ and /j+i crosses Xj.iXj. In 
this case, however, each such pair of neighbours is contained in the trap region generated 
by any other one and hence we get a cycle of order 4 consisting of some Xj, yi-i+i, Xj, yjj+i, 
where is the intersection point of /j and /j+i and likewise for yjj+i- In either situation 
Q~^[li,xi; ... ; /fe, Xk] contains a cycle as required. This completes the proof of Theorem 5. □ 
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5.3 Full characterisation of collective factorisation for rectangu- 
lar fields 

As we already have mentioned above, in the particular case of rectangular Markov fields our 
knowledge is more complete than in the general setting. In fact, we are able to show that 
the necessary precedence graph acyclity condition for collective factorisation on all scales, 
as stated in combined Lemmas 1 and 2, becomes in this context a sufficient condition as 
well. This is made precise in the following theorem. 

Theorem 6 Assume that the field A'^ is rectangular and let {l-hXi)^^^^ Xi G /j, be a 
collection of pairwise distinct lines and points in general position. Then the following are 
equivalent: 

1. The precedence graph Q[li,Xi] . . . ; /fc, x^] is acyclic. 

2. The collection {li,Xi) admits collective factorisation of correlations on all scales. 

Proof The implication from 2. to 1. follows directly from Lemmas 1 and 2. Thus, it only 
remains to establish the implication from 1. to 2. To this end. we make first the following 
observation, which is also of its own intrinsic interest in the context of rectangular fields. 

Given two non-parallel lines I and on whose parallel translates the activity measure 
M. is concentrated, we say that a set D C is {I, Z'}-convex if for any two points x,y E D 
lying on a common translate of I or I' the entire segment xy is contained in D. Clearly, {I, l'}- 
convexity is a much weaker concept than the usual convexity, in particular a {I, r}-convex 
set does not even have to be connected. We shall write conv {i^i'}{A) for the {Z, Z'}-convex 
hull of A C R^, that is to say the smallest {I, r}-convex set containing A. Now, the crucial 
observation is that in the context of the generalised dynamic representation in Section 4 
specialised for the considered rectangular field, the convexity assumption imposed on the 
growing family (Dt) there can be relaxed to the {I, r}-convexity here without any further 
modifications of the theory. Moreover, instead of having the initial set Dq consist of a single 
point, we can now let it consist of any positive finite number of points under the {/, /'}- 
convexity requirement. Indeed, the proof of Theorem 3 carries over verbatim under these 
relaxed conditions on {Df) and the {/, /'}-convexity of {D^) ensures that the intersections of 
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the field lines with (Dt) are always connected. It is essential to note at this point that the 
convexity requirement cannot be lifted for the target domain D though. The point is that 
the consistency statement in Theorem 2 does essentially require convexity of the domain 
of the field in its proof and, consequently, a rectangular field constructed in some non- 
convex but {/, /'}-convex set may fail to coincide with the corresponding restriction of the 
whole-plane consistent field. This phenomenon is closely related to the rather unintuitive 
fact that there exist pairs of compacts Di C D2 ^ consisting of a {/, /'}-convex set 
Di and a convex set D2 with the property that Di cannot be extended to D2 through 
a Hausdorff-continuous growing family of {/, /'}-convex compacts, in sharp contrast to 
the case of convex sets where such an extension is always possible. These issues are not 
discussed in further detail here as falling beyond the scope of the present article, yet they 
are a subject of our active research in progress because we believe that the knowledge of 
the geometry of such non- extensible sets and their corresponding defective fields may shed 
some further light upon the higher-order correlation structure of A''^. 

To proceed, recall that the vertex set of the precedence graph Q[li,xi; . . . ;lk, x^] con- 
sists of the generating points Xj's and of Ui/s arising as the intersection points of the 
corresponding li and Ij. Since the precedence graph is acyclic, there exists a complete or- 
dering of its vertices compatible with the structural partial order induced by directions of its 
edges. In addition, in view of the minimality of Xj's in the structural order, we can assume 
that no y.^.-vertex (intersection vertex) precedes an x.-vertex (generating vertex) in the con- 
sidered complete ordering, that is to say the ordered collection of vertices can be written as 
{xi, . . . , Xk', vi, . . . , Vm) where vi, . . . ,Vm are all y.,. -vertices of Q[li,xi; ... ;lk, Xk\- Say that 
a vertex u of the precedence graph is a direct predecessor of another vertex w if there is a 
directed edge from u to w containing no other vertices of the graph, which amounts to di- 
rect precedence of u over w in the structural order on the vertices of the precedence graph. 
Note that y.,. -vertices have precisely two direct predecessors each, whereas x. -vertices have 
no predecessors. We construct inductively a growing family {Dt) of /'}-convex sets in 

by putting Dt := convj;^;/} D^, where 

• We put Dq = {xi,... , Xk}. 

• For t — i/m, 1 < i < m we set D° = {xi, . . . , x^, Vi, . . . , Vi}. 

• For t e {{i — l)/m, i/m), 1 < i < m we let wi,W2 be the direct predecessors of the 
vertex Vi (note that tfi, £ for i < (i — l)/m by the construction) and we define 
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as the union of 

— 77° 

— the point wi + m(i — {i — l)/m){vi — wi), 

— the point W2 + m{t — {i — l)/m){vi — W2). 

In other words, initially the set Dq consists of {xi, . . . , Xk}- Since the collection (Zj, Xi) is in 
general position, we clearly have Dq — conv^i i/y Dq — conv{;_;/}{xi, . . . , Xk} — {xi, . . . , Xk}- 
Next, the intersection points Vi are subsequently added to the set D^, in time intervals of 
length 1/m. With the intersection points added, the {Z, r}-convex hull of the increasing 
vertex collection does no more coincide with the collection itself. The time between the 
moments of y.,. -vertex additions is therefore used to interpolate the family (D^), thus 
keeping the growth of its {I, r}-convex hull (Df) Hausdorff continuous. Indeed, the growth 
of {D^) is Hausdorff continuous by definition and it is easily seen that each new point in 
Dt+dt \ Dt arises as a shift of a point in D^_^_^^ \ along a translate of either / or I' and, 
consequently, Dt+dt \Dt is a, union of two rectangles of infinitesimal width 0{dt) built on 
two possibly degenerated segments parallel to / and respectively. Note that in general 
this would not be the case if fj's were not ordered compatibly with the structural ordering 
- adding a single new intersection point to might then possibly result in some further 
intersection points, not yet present in D^, falling into its {/, /'}-convex hull Dt, which might 
in its turn produce an extra non-degenarate rectangle in Dt thus violating the Hausdorff 
continuity. The importance of imposing on f j's an ordering compatible with the structural 
order for (Zj, Xi)^^^ is that it ensures that the intersection points f j in Dt are precisely those 
present in D^ with no extras, that is to say for t e [i/m, (i + l)/m) fl [0, 1], i = 0, 1, . . . , m, 
we have 

{Xi, ... ,Xk,Vi,... , Vm} r]Dt = {Xl, ... ,Xk,Vi,... ,Vi}nD° = {Xi,... , Xk, Vi,... , Vi}. 

In particular, the time order in which the points of the precedence graph Q[li,xi; . . . ;lk,Xk\ 
show up in Df is compatible with the structural ordering for the collection {li,Xi)f^i. 
Observe as well that Di is a convex set (a rectangle in fact) containing all Xj's for i — 
1, . . . ,k. Moreover, even though (Dt) does not satisfy the condition (D4) in Section 4, 
it can be easily modified by local smoothing to yield an increasing family (Dt) of {/, /'}- 
convex sets satisfying the conditions (D1),(D3),(D4),(D5) in Section 4 and enjoying the 
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property (inherited from (Dt)) that A(Zj) = Xi where A(-) is the anchor mapping induced 
by Df. Furthermore, by the same local smoothing modification Di can be taken convex 
and containing all Xj's in their interior. In view of the discussion on the particular form 
of the generalised graphical construction for rectangular fields and taking into account the 
(GErLineBirth) rule of this construction we conclude the required collective factorisation 
of correlations on all scales, thus completing the proof of Theorem 6. □ 

It should be emphasised at this point that the above argument apparently cannot be 
repeated for more general fields. Indeed, with all directions allowed the condition (D2) 
cannot be lifted and Dq has to be a singleton which makes our construction above break 
down. Moreover, in case when only a finite number of, but more than two, directions are 
allowed, the construction breaks down as well - even though a discussion parallel to the 
above can be provided based on the concept of convexity in the directions of the field, 
which does again allow to relax the offending condition (D2), wc may in general loose the 
Hausdorff continuity of Dt at the moments when two parts of Dt growing in two different 
directions become for the first time connectablc along some other direction of the field 
(never to happen when there are only two directions allowed). 

5.4 Integral correlation measures 

With the knowledge of the correlation structure of the field as provided so far in this 
section, we are now in a position to determine natural integral correlation measures for 
the field. To this end, for a bounded open convex D CM? and for a planar graph in D (a 
collection of edges, not necessarily in Fd) we construct its first order directional measure 

T>^'\l]--= E ^(^)^ne] (13) 

eGEdges(7) 

with £{■) denoting the usual Euclidean length and with 5;[e] standing for the unit mass at 
l[e] e [[D]]. The so-defined V^'^\-f) is a finite purely atomic measure on [[D]]. Likewise, we 
consider the second order directional measure on [[D]] x [[D]] given by 

^^'^[7]:= E E ^(ei)^(e2)5(,[e,],ne.]). (14) 

ei6Edges(7) 62 6 Edges (7) 

We shall use the expectations of these measures with 7 drawn from as natural char- 
acteristics of the directional nature of the field. To this end, we put 

:=EPW[>1^]. (15) 
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Write also 

:= EP(2) [^^] -A^®A^ (16) 

where (g) denotes the measure producting operation. Note that E;^ has its natural inter- 
pretation as the directional covariance measure. Below we explicitly determine both the 
first and second order directional measures for A'^. 

CorollEiry 1 For I e [[D]\ we have 

A^{dl) =£{lr\D)M{dl). 

Proof For each / e [[D]] we have 



A-^{dl)= / a'^[dl,x]i{dx) 
JlnD 

where £ is the length element. Our statement follows now by Theorem 4. □ 

Corollary 2 For I1J2 G [[D]] we have 

Y^M^i dl ) = ! IinDlinD^M-'^Mi[[xy]]))e{dx)e{dy)M{dl), if h ^k^l, 
^ ^ ^' \ 0, otherwise. 

Proof For /i,/2 G [[D]] we have 

T.^{dh,dh)= I I a^[dh,Xi-dl2,X2Y{dx2)l{dxi)- A^{dh)A^{dh). 
JhnD JhnD 

Our assertion follows now by Theorem 4. □ 

5.5 Alternative correlation measures 

In addition to the events that there exist field lines passing through given points in given 
directions as considered in the definition of . .] above, one can also study events that 
there exist field angles with vertices at given points and with their arms extending along 
given half-lines. In formal terms, we put 

^^[dZi; . . . ; dZk] := dF {^l^Z, e Angles(^^)) , (17) 

where Angles(7), 7 G F, denotes the collection of convex angles in the admissible polygonal 
graph 7 and where Zj := Z[li, I'i], i — 1, . . . , k, stands for the convex angle between two 
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directed half-lines k and I'i outgoing from the vertex of the angle and directed away from 
the vertex. It should be noted at this point for formal completeness that the concave angles 
arising as the complements of the convex angles of 7 are not included into Angles (7). In 
analogy to Lemma 3, a similar factorisation result for the probabilities of such events can 
be concluded from the generalised dynamic representation. 

Corollary 3 Assume the collection Zi, . . . , Zfc is such that the vertices Xi of Zi are all 
different and such that the interior of the convex angle Zi is disjoint with the convex hull 
conv({a;2, . . . , Xk})- Then 

<;^[dZ,; . . . ; dZ,] = M{dh)M{dl[)<;^[dZ2; . . . ; dZ,]. 

Proof Indeed, under the assumptions of the corollary the increasing family Df in the 
generalised dynamic representation in Section 4 can be chosen so that it first hits all 
X2, ■ ■ ■ ,Xk and finally it reaches Xi before hitting any other point of Zi. Our claim follows 
now by the rule (GErVertexBirth) of the generahsed representation. □ 

Suppose now that we assign alternating labels +1 and —1 to regions separated by the 
contours of the field which thus become interfaces between +1 and —1 phases. This gives 
rise to two possible label assignments in IR^ and we pick any of them with probability 1/2. 
Write [^-^ja; for the label assigned by at a; e M?. Then for two points |/ G we 
can define the two-point label-correlation function 

P^y m^]M^]y - m'^im'^]y (18) 

It is easily seen that, writing n^^y '■— caxd{xy O A-^) for the number of edges of the field 
A^ crossing the segment xy, 

P^y = ^{^x,y is even ) - "^{n^^y is odd ). 

RecaUing from the Lineeir sections statement of Theorem 2 that n^^y is Poisson with 
parameter with xy denoting the segment joining x to y, we conclude that 

CoroIIairy 4 We have 

P^y = eM-Mm]]))- 
Roughly speaking, this shows that for regularly behaved activity measures M. the two- 
point label-correlation functions of the field A'^ exhibit exponential decay with the distance 
between x and y which places the inverse temperature /3 = 1 corresponding to the consistent 
regime in the high temperature regime for polygonal fields, above the phase transition point. 
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6 Disagreement loop death and birth dynamics 



In this section we discuss a random dynamics on the space Fd of admissible polygonal 
configurations which leaves the law of the consistent field invariant, with D standing 
as usual for an open bounded convex set in R^. We also provide a modification of this 
dynamics which can be used for Monte-Carlo simulation of A^'^ for all /3 e IR. We build 
upon [15] in our presentation of the standard dynamics based on an important concept of 
a disagreement loop. 

6.1 Standard disagreement loop dynamics for consistent regime 

To proceed we place ourselves within the context of the standard dynamic representation 
as given in Subsection 3.1 and suppose that wc observe a particular realisation 7 G F^i 
of the polygonal field and that we modify the configuration by adding an extra birth 
site xq to the existing collection of birth sites for 7, while keeping the evolution rules 
(El-4) for all the particles, including the two newly added ones if Xq ^ D and the single 
newly added one if xq G dD. Denote the resulting new (random) polygonal configuration 
by 7 © Xq. A simple yet crucial observation is that for xq E D the symmetric difference 
7A[7©a;o] is almost surely a single loop (a closed polygonal curve), possibly self-intersecting 
and possibly chopped off by the boundary. Indeed, this is seen as follows. The leftmost 
point of the loop 7 A [7 ® xq] is of course xq. Each of the two new particles pi,p2 emitted 
from xq move independently, according to (El — 4), each giving rise to a disagreement 
path. The initial segments of such a disagreement path correspond to the movement of 
a particle, say pi, before its annihilation in the first collision. If this is a collision with 
the boundary, the disagreement path gets chopped off and terminates there. If this is 
a collision with a segment of the original configuration 7 corresponding to a certain old 
particle pa, the new particle pi dies but the disagreement path continues along the part of 
the trajectory of which is contained in 7 but not in 7 ® xq. At some further moment 
P3 dies itself in 7, touching the boundary or killing another particle p4 in 7. In the second 
case, however, this collision only happens for 7 and not for 7©xo so the particle p4 survives 
(for some time) in 7 © xq yielding a further connected portion of the disagreement path 
initiated by pi, which is contained in 7 © xq but not in 7 etc. A recursive continuation of 
this construction shows that the disagreement path initiated by pi consists alternately of 
connected polygonal subpaths contained in [7 © xq] \ 7 (call these creation phase subpaths) 
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and in 7 \ [7 ® xq] (call these annihilation phase subpaths) . Note that this disagreement 
path is self-avoiding and, in fact, it can be represented as the graph of some piecewise 
linear function t 1— > y{t). Clearly, the same applies for the disagreement path initiated 
by p2- An important observation is that whenever two creation phase or two annihilation 
phase subpaths of the two disagreement paths hit each other, both disagreement paths die 
at this point and the disagreement loop closes (as opposed to intersections of segments of 
different phases which do not have this effect). Obviously, if the disagreement loop does 
not close in the above way, it gets eventually chopped off by the boundary. We shall write 
A®[a;o;7] = 7^[7 © ^0] to denote the (random) disagreement loop constructed above. It 
remains to consider the case xq G dD, which is much simpler because there is only one 
particle emitted and so A®[a;o;7] = 7A[7 © Xq] is a single self-avoiding polygonal path 
eventually chopped off by the boundary. We shall often abuse the language calling such a 
disagreement path A®[a;o;7] a (degenerate) disagreement loop as well. 

Likewise, a disagreement loop (or path) arises if we remove one birth site Xo from the 
collection of birth sites of an admissible polygonal configuration 7 G F^i, while keeping 
the evolution rules for all the remaining particles. We write 7 for the configuration 
obtained from 7 by removing xq from the list of the birth sites, while the resulting random 
disagreement loop (or path) is denoted by A®[a;o;7] so that A®[xo;7] — 7A[7Qa;o]. We 
refer the reader to Section 2.1 in [15] for further discussion. 

With the above terminology we are in a position to describe a random dynamics on 
the configuration space To, which leaves invariant the law of the consistent polygonal 
process A^. Particular care is needed, however, to distinguish between the notion of time 
considered in the dynamic representation of the field as well as throughout the construction 
of the disagreement loops above, and the notion of time to be introduced for the random 
dynamics on F^) constructed below. To make this distinction clear we shall refer to the 
former as to the representation time (r-time for short) and shall reserve for it the notation t, 
while the latter will be called the simulation time (s-time for short) and will be consequently 
denoted by s in the sequel. 

Consider the following pure jump birth and death type Markovian dynamics on To, 
with 7s = 7^* standing for the current configuration 

DL:birth With intensity {{M)){dx)ds for x G -D and with intensity {{M;dD)){dx)ds 
for X e dD set ^s+ds 7s ® x, 
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DLrdeath For each birth site x in 7s with intensity ds set ^s+ds '■— Is © x. 

If none of the above updates occurs we keep ^is+ds = 7s- It is convenient to perceive the 
above dynamics in terms of generating random disagreement loops A and setting 7^+^^ := 
7sAA, with the loops of the type A®[-, •] corresponding to the rule DL:birth and A®[-, •] 
to the rule DL:death. 

As a direct consequence of the dynamic representation of the consistent field A-^ we 
obtain 

Theorem 7 The distribution of the polygonal field is the unique invariant law of the 
dynamics given by DLrbirth and DLrdeath. The resulting s-time stationary process is 
reversible. Moreover, for any initial distribution of 70 the laws of the polygonal fields 7^ 
converge in variational distance to the law of A^ as s ^ 00. 

The uniqueness and convergence statements in the above theorem require a short justifi- 
cation. They both follow by the observation that, in finite volume, regardless of the initial 
state, the process 7^ spends a non-null fraction of time in the empty state (no polygonal 
contours). Indeed, this observation allows us to conclude the required uniqueness and con- 
vergence by a standard coupling argument, e.g. along the lines of the proof of Theorem 
1.2 in [13]. 

6.2 Standard disagreement loop dynamics for general tempera- 
tures 

Below, we show that the laws of the Gibbs-modified polygonal fields A'^''^, G M, arise as 
the unique invariant distributions for appropriate modifications of the reference dynamics 
DLrbirth, DLrdeath. The main change is that the birth and death updates are no 

more performed unconditionally, they pass an acceptance test instead and are accepted 
with certain state-dependent probabilities whereas upon failure of the acceptance test the 
update is discarded. Consider the following dynamics 

DL : birth[/3] With intensity {{M)){dx)ds for x e D and with intensity {{M; dD)){dx)ds 
for X e dD do 

• put 5 := 7s © X, 

• accept 6 with probability exp(-[/3 - 1]{L^{S) - L^(7s))) if L^{S) > L^{%) 
and with probability 1 otherwise. 
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• if accepted, set js+ds '■— otherwise keep 7s+ds '■— 7s- 
DL : death[/?] For each birth site x in 7^ with intensity ds do 

• put 5 := 7s X, 

• accept 6 with probabihty exp{-[(3 - 1]{L-^{S) - L^i'ys))) if L^{S) > L^{ls) 
and with probabihty 1 otherwise, 

• if accepted, set ^s+ds '■— otherwise keep 7s+ds '■— ^s- 

In analogy with its original reference form DL:birth, DL:death, the above dynamics 
should be thought of as generating random disagreement loops A and setting ^s+ds '■ = 
7AA provided A passes the acceptance test. The following theorem justifies the above 
construction. 

Theorem 8 For each P E'^ the law of the polygonal process A^'^ is the unique invariant 
distribution of the dynamics DL : birth[/3], DL : death[/9]. The resulting s-time stationary 
process is reversible. For any initial distribution of 70 the laws of the polygonal fields 7^ 
converge in variational distance to the law of A^'^ as s ^ 00. 

Indeed, the invariance follows by a straightforward check of usual detailed balance con- 
ditions whereas the convergence and uniqueness statements are verified as in Theorem 
7. 

6.3 Generalised disagreement loop dynamics 

Even though the above disagreement loop dynamics in its standard version has been con- 
structed using the standard dynamic representation of Subsection 3.1, exactly the same 
can be made for the generalised dynamic representation of Section 4. As easily verified, 
this also leads to disagreement loop and path creation and annihilation, the only difference 
being that boundary birth events in the standard set-up are replaced by general line birth 
events in context of the generalised dynamic construction. Vertex creation and annihilation 
yield disagreement loops whereas line creation and annihilation yield disagreement paths. 
Thus, the DL: birth and DL: death moves get replaced by 

GenDL:birth With intensity {{Ai)){dx)ds foi x E D set ^s+ds '■= 7s ©a;. With intensity 
M.{dl)ds for / G [[D]] set js+ds '■= 7s © ^ where 7^ © / arises from 7s by letting / 
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be born at its anchor point A{1) and making it extend thereupon according to the 
evolution rules of the generalised construction, which gives rise to a disagreement 
path. 

GenDL:death For each point birth site x in 7^ with intensity ds set ^s+ds '■— Is Qx. For 
each line birth site of a line Hn 7^, with intensity ds set 7s+ds '■— IsQl-i where 
again 7^ G I arises from 7^ by kiUing at A(Z) the field line I of 7^, which gives rise to 
a disagreement path. 

In full analogy with Theorem 7 we have 

Theorem 9 The distribution of the polygonal field A'^ is the unique invariant law of the 
dynamics given by GenDLibirth and GenDLrdeath. The resulting s-time stationary 
process is reversible. Moreover, for any initial distribution of 70 the laws of the polygonal 
fields 75 converge in variational distance to the law of as s ^ 00. 

The modification GenDL[/3] of the GenDL dynamics for general /? e M goes along exactly 
the same fines as the construction of DL[/9] from DL, by introducing suitable acceptance 
tests, and hence we omit the standard details here. Obviously, the corresponding version 
of Theorem 8 holds in full analogy. Clearly, the so-defined dynamics depends then on 
the choice of the increasing family D^. Note that hybrid dynamics can also be considered, 
switching between different choices of Dt in the course of the simulation time s, which is 
sometimes beneficial e.g. in image processing applications, see [17]. 

6.4 Defective disagreement loop dynamics 

To conclude the discussion of disagreement loop dynamics we introduce one more variant 
thereof, which we call defective dynamics. The idea stems from the observation that the 
standard dynamic construction in Subsection 3.1 and its generalised version in Section 4 
rely respectively on the time axis direction and on the anchor mapping A : [[D]] D 
to decide in which of the two possible directions to proceed along a new line each time a 
directional update occurs and a new edge is to be created. This decision mechanism is then 
inherited by the disagreement loop dynamics where the current direction of an unfolding 
disagreement loop changes in the following situations: 

1. upon a directional update during a creation phase. 
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2. when a directional update is present on the path being erased during an annihilation 
phase, 

3. when switching between creation and annihilation phases. 

In all these situations the choice of the direction to follow along the new hne is always 
the same: forward in time in the standard case and away from the corresponding anchor 
point in the generalised case. In certain contexts there may be good reasons to try to 
overcome these restrictions and to allow for more general directional decision strategies 
even at the cost of obtaining disagreement loop dynamics which do not correspond to any 
well-defined dynamic representations. One of such situations arises when studying the 
correlation functions of the polygonal field as defined in Section 5 and we shall indeed use 
the defective construction to provide below a proof of exponential correlation decay for 
rectangular fields, see Section 7. Another setting where more general directional decision 
strategies may be of use are image processing applications of polygonal fields, see [17]. 

The idea underlying the defective disagreement loop dynamics is hence to keep the cru- 
cial feature of the above DL and GenDL dynamics, being that the symmetric difference 
between a configuration and its update is always a disagreement loop algorithmically con- 
structible in the course of alternating creation and annihilation phases, but to give up the 
requirement that the disagreement loop generating mechanism be determined by a dynamic 
representation. Such an approach makes this disagreement loop generating mechanism the 
core of the construction. In our construction of the defective dynamics we shall assume 
that upon each directional update along an unfolding disagreement loop, the decision on 
which of the two possible directions to choose along the new line is taken as a determin- 
istic function of the current portion of the disagreement subpath built by the time of the 
update. To put it in formal terms, assume we are given a measurable directional decision 
mapping/rule 4) — 4>d assigning to each connected non-closed but possibly self- intersecting 
polygonal path A contained in the domain D of the field and to a line I e [[D]] passing 
through one of the two ends of A, one of the two possible directions along I. In our dynamics 
below (f>{X,l) will be used to determine the direction in which the disagreement path A is to 
extend upon a directional update into /. For a polygonal configuration 7 G F/j, a straight 
line / G [[D]] assumed not to be an extension of an edge of 7 and a point a: G / \ 7 we 
construct the configuration 7 © (/, x) G F/j, interpreted as the result of adding to 7 a linear 
germ (/, x) born at x and initially evolving along I 3 x oi, in other words, creating a new 
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field fine for 7 at x. The construction goes as follows. We add a fine birth site for I at x 
at the time t — whereupon we let the newly born linear germ extend in both directions 
along I away from x, say with constant unit speed for definiteness although other options 
are also possible as discussed below. While extending, the resulting new edge is subject 
to directional updates: for each line I' which the growing edge crosses the probability that 
a directional update occurs and a new growth direction along /' is assumed is M{dl') - if 
this happens at some time moment the new direction along I' is always chosen according 
to the directional decision rule 0(A(, /') and the growth continues with constant unit speed, 
with A( standing for the current path constructed by the time t. At some point, the result- 
ing polygonal path may hit the boundary of D or an existing edge e of 7, in which case 
it stops growing in this direction. Upon a collision with an edge we switch from creation 
mode to annihilation mode and proceed (say, again with constant unit speed) from the 
collision point along the edge e in the direction indicated by 0(At,/[e]), erasing this part 
of the edge. This continues till we reach the boundary dD or another edge e' of 7. In the 
latter case, if the edge e' extends in the direction agreeing with 0(At, l[e']) from the collision 
point, we proceed along e' and keep erasing it. If, however, the edge e' extends from the 
collision point in the direction opposite to 0(At, /[e']), we switch back to the creation mode 
and start extending e' in the direction indicated by 0(At, l[e']). We continue this procedure 
for both branches of the linear germ emitted from the line birth site x). Both branches 
give rise to directed polygonal disagreement paths constituting the symmetric difference 
7A(7 ® {l,x)), both consisting of alternating polygonal subpaths corresponding to cre- 
ation and annihilation stages of the construction and respectively referred to as creation 
phase and annihilation phase subpaths in the sequel. The procedure terminates when both 
branches get killed in collision with dD. It may be useful to observe that the requirement 
that both branches of the disagreement loop grow with constant unit speed has no ef- 
fect on the final shape of individual branches but it determines the precedence of their 
intersection points, should they meet. Thus, in principle, the unit growth speed may be 
alternatively replaced by other natural options such as letting the right (left) branch grow 
first till it reaches the boundary and to let the left (right) branch grow thereupon, or even 
admitting different growth speeds for both branches or speeds depending on the current 
form of the disagreement loop under construction. Below we shall always assume one of 
these options is chosen and fixed, under each such choice the discussion below remains 
valid. The co-existence of two disagreement branches obliges us to impose one more rule 



36 



on the construction of 7 ® (Z, x) : if a creation phase subpath of one branch hits a creation 
phase subpath of another branch, both branches get killed at this point, referred to as the 
cut-off point below, their continuation possibly already constructed gets chopped off and 
we consider the disagreement loop complete thus terminating the construction; likewise for 
two annihilation phase subpaths of different branches colliding. If two subpaths of different 
phases coming from different branches meet, nothing happens though. Note that there can 
be several such cut-off points and the disagreement loop construction terminates at the 
first cut-off point with respect to the ordering imposed by the time flow - consequently 
the choice of a particular growth speed protocol as discussed above may affect the choice 
of the loop-closing cut-off point and hence also the final shape of the disagreement loop. 

So far, this construction looks very similar to the construction of the disagreement 
paths and loops in DL and GenDL above. There is one crucial difference though: it 
may happen here that a self-intersection occurs along one branch of the disagreement path 
being constructed - this could never have happened in the disagreement loop dynamics DL 
and GenDL corresponding to dynamic representations because of the temporal ordering 
along branches of disagreement paths/loops there, and this is the price we pay here for 
the freedom of choice of the directional decision mapping ■). If such a self-intersection 
occurs past the cut-off point, we just ignore it. If it happens before the cut-off, we are in 
trouble though. A look at the description of the disagreement path construction procedure 
here readily shows that the presence of such a self-intersection and hence of a cycle along 
a single disagreement branch before the cut-off means 

• either that at this point we start discarding the changes previously introduced by 
erasing segments previously created (meeting of two creation phases) 

• or that the ®{l,x) operation is not reversible (creation phase edge hits annihilation 
phase subpath preceding it along the disagreement branch) in that the disagreement 
loop 7A(7 ® (Z, x)), even though it can be used for the transition 7 — > (7 ® {I, x)), 
is no more a valid disagreement loop upon switching the creation and annihilation 
phases and thus cannot be used for the reverse transition (7 © (/, x)) — > 7. 

The remaining two situations: intersection of two annihilation subpaths or an annihilation 
edge hitting a creation subpath preceding it along the disagreement branch, cannot occur 
in our construction as this would either imply the existence of a node of order higher 
than 2 in the original configuration 7 (in the first case) or stand in contradiction to the 
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disagreement loop construction rules (in the second case) . To prevent the above pathologies 
from happening we say that the construction of the disagreement loop/path fails if a cycle 
occurs along one of the branches. Consequently, the configuration 7 © {l,x) is not always 
defined, hence the name defective dynamics. 

In full analogy with the above definition of the ® operation, we can define the configu- 
ration 7 (/, x) e To for line / extending an edge e e Edges(7), 7 e F/), with x e e, which 
is to be interpreted as initiating the annihilation phase in 7 from the linear germ (/, x) or, 
in other words, killing/annihilating the field line / of 7 at the point x. This operation also 
gives rise to a disagreement loop and its outcome may be undefined as well, on equal rights 
with 0. 

The operations 0(/, x) and 0(/, x) combine the properties of line birth/death and vertex 
birth/death operations of the generalised disagreement loop dynamics in Subsection 6.3 
and there are no separate 0/, 0a;, 0/, 0a; operations in the defective dynamics. 

Directly from the definition of the line creation operation it follows that if for some 
7 e Vjy and a linear germ {l,x), a; G / \ 7, another configuration 7' e Fj) is reachable from 
7 upon adding the hnear germ (/, x), that is to say 7' is a possible outcome for 7 (/, x), 
then the probability element that 7 ® (i, x) = 7' is 

exp(-L-^(7'\7)) n M{dl[e]), (19) 

eGEdges(7'\7), e^x 

where Edges(7' \ 7) stands for the collection of edges of 7' which are not extensions of 
edges of 7 and thus they have to result from directional updates during creation phase. 
Indeed, the product ]^ (rf/fe]) in (19) comes from the probability cost of creation phase 
directional updates yielding all edges in Edgcs(7'\7) save the initial one germinating from 
X, whereas exp(— L-^(7'\7)) is due to the absence of directional updates along the created 
polygonal paths of 7' \ 7. In a similar way we see that if for some {I, x) the configuration 
7 e Fi3 contains an edge e with x & e and I — l[e] then for another configuration 7' 
reachable from 7 by annihilating the line Z at x we have the probability element that 
7 G {I, x) — 7' given by 

exp(-L^(7'\7)) n M{dl[e]). (20) 

eeEdges(7'\7) 

It is also crucial to observe that each valid transition 7 ^ 7' = 7 (/, x) admits its inverse 
obtainable as 7' — > 7 = 7'© {I, x) under the same directional decisions made in the course of 
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the disagreement loop construction but with the creation and annihilation phases switched 
- this is due to the fact that the directional decision rule •) only takes into account the 
current state of the unfolding disagreement loop (which is the same for 7 — > 7' and 7' ^ 7 
as the symmetric difference) rather than the actual configuration (which differs for both 
transitions) . 

Consider now a measurable set C [[D]] and a measurable mapping A* : — > D 
such that A*(/) e I. These are interpreted respectively as the set of lines allowed for cre- 
ation/annihilation and the assignment to each line of its germination point. The mapping 
A* can also be randomised as long as it is kept independent of the dynamics below, but 
for simplicity we do not discuss the details of this straightforward extension as it is not 
going to be used in this paper. Define the following defective disagreement loop dynamics 
unfolding in s-time: 

DefDLrbirth With intensity M{dl)ds for I e Al set -fs+ds := 7s © {l,A*{l)). Should the 
proposed update fail, keep js+ds — Is- 

DefDL:death For each line I G Al extending an edge e of 7^ such that A*(/) G e, 
with intensity ds set js+ds '■— 7s © (^)-^*(0)- Should the proposed update fail, keep 

7s+(is 7s- 

Recalling the above transition probabibilities (19) and (20) corresponding to ® and 
moves and taking into account the discussion following their displays allows for a straight- 
forward verification of the detailed balance conditions which readily yields 

Theorem 10 The distribution of the polygonal field A"^ is an invariant law for the dy- 
namics given by DefDLrbirth and DefDLrdeath. The resulting s-time stationary process 
is reversible. 

Note that in strong contrast to analogous Theorems 7 and 9 wc do not state uniqueness or 
convergence to the law of the polygonal field A^^ here and this is because in general neither 
of them holds. Indeed, this is due to the fact that a defective dynamics does not have to 
be transitive: there may exists pairs of admissible polygonal configurations which cannot 
be reached from each other by © and moves. This is closely related to the possible 
failures of these updates and it is the price to pay for the free choice of the directional 
decision mapping. The structure of the recurrent classes of defective dynamics on 
is a subject of our research in progress. Observe as well that, in view of the discussion 
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following the displays (19) and (20), the requirement that the directional decision rule 
0(-, •) only depends on the current state of the disagreement loop rather than on the entire 
actual configuration seems indispensable for the detailed balance and the reversibility of 
the dynamics. 

Veiriants and applications of the defective dynamics The modification DefDL[/3] 
of the defective disagreement loop dynamics DefDL for general /3 7^ 1 is straightforward 
and goes by introducing suitable acceptance test for update proposals in full analogy with 
the corresponding modification of the DL dynamics aimed at obtaining DL[/3] in Subsec- 
tion 6.2. 

The class of possible directional decision rules is extremely rich and various strategies 
can be chosen encompassing different goals. One possible option of constructing the di- 
rectional decision mapping goes by choosing an arbitrary measurable generalised anchor 
mapping A : [[D]] D with A(/) G / to be interpreted as the initial point of the line /. 
Then the direction 0(A, /) the disagreement loop A is to assume for its unfolding along / 
upon a directional update is always set away from the anchor point A(/). Note that if the 
anchor mapping corresponds to an increasing family (D^) of convex compacts satisfying 
(Dl-5) as in Section 4 then the resulting defective dynamics © and operations coincide 
respectively with the line birth and line death operations of the generalised diasgreement 
loop dynamics described in Subsection 6.3. 

The directional decision mapping can also be randomised, as long as 0(At, •) is indepen- 
dent of the current field configuration conditionally given A^. Such randomised strategies 
will not be used in this paper though. 

Due to its flexibility, we envision manifold applications of the defective dynamics in 
our research in progress on polygonal fields. In the first place we are planning to use it to 
gain further knowledge on the higher order correlation structure of consistent polygonal 
fields. To illustrate the techniques we are currently developing we shall employ the defective 
dynamics to establish exponential decay of correlations in the particular case of rectangular 
fields in Section 7 below. Another important field in which we anticipate the use of this 
dynamics is digital image segmentation - the defective dynamics will provide a rich and 
fiexible class of new Monte-Carlo moves for polygonal field-based Bayesian segmentation 
algorithms we developed in joint work with Kluszczyhski and van Lieshout [10, 11, 12, 17]. 
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7 Exponential decay of dependencies for consistent 
rectangular fields 



The purpose of this section is to make use of the defective disagreement loop dynamics 
developed in Subsection 6.4 above in order to establish exponential decay of higher order 
correlations for the consistent rectangular field A'^. To simphfy the presentation we assume 
throughout this section with no loss of generality that the admissible directions 1,1' of the 
field are parallel to the coordinate axes. Assume also that a collection (Zj, Xi)^^i, Xi e li, in 
general position as in Section 5 is given such that xi is its topmost point with the highest 
y coordinate. Further, let TZ[x2, . . . ,Xk] he the smallest rectangle with sides parallel to 
coordinate axes containing all remaining points X2, ■ ■ ■ ,Xk of the considered collection. 
We know from Lemma 3 that if does not hit the convex hull of X2, ■ ■ ■ ,Xk then the 
correlation function a-^[li,xi; l2,X2; ■ ■ ■ ;lk, Xk] factorises as o"-^[/i, Xi]a^[l2, X2\ ■ ■ ■ 'ilk-, Xk\- 
Here in Lemma 4 we show that also when li does hit conv({a;2, . . . iXk}) and hence it 
also hits the rectangle TZ[x2-, ■ ■ ■ , a^^j], the correlation Xi\l2,X2\ . . . ;lk, Xk] approaches 

a-^[li, Xi]a-^[l2, X2; ■ . ■ ',lk,Xk] exponentially fast in the distance measure between xi and 
TZ[x2, ... ,Xk] along h. 

Lemma 4 With notation as above, if h hits TZ[x2, ... ,Xk\ then we have 
a^lh, xi, I2, X2;...; k, Xk] = cr^ih, Xija^lk, X2, . . . ; Xk] 

where li[xi ^ 1^1x2, ... ,Xk]] is the segment of li between Xi and TZ[x2, ... ,Xk\- 

Proof Our proof relies on estimating the conditional probability of the event 

S[dli,xi] := {3eeEdges(.A-^) xi e e, l[e] G dh} 

given the event 

S[dl2, X2]... ; dlk, Xk] := {Vf=23eeEdges(.4^) Xi e e, l[e] e dk}. 
We shall show that 

P {S[dh, xi] \S[dl2, X2;...; dlk, Xk] ) = M{dh) (1 + (e-Millhl-i^n-2,...,x,]]]])^^^ ^ 

(21) 
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which will complete the proof in view of (10) [recall that Xi e k, i — 1, . . . ,k] and Theorem 
4. To proceed, denote by a rectangular field coinciding in law with conditioned 
on the event S[dl2, X2; ■ ■ ■ ;dlk,Xk\- Note that since li crosses TZ[x2, ... ,Xk\, it has to be 
a vertical hne. Construct the following directional decision rule ^(•, •) for the defective 
dynamics: 

• Whenever turning into a horisontal direction, move away from li, that is to say 
rightwards if the turning point lies to the right from li and leftwards if the turning 
point lies to the left of li. When the turning point lies on li, move rightwards. 

• Whenever turning into a vertical direction, move upwards. 

Moreover, let D :— TZ[xi, X2, ■ ■ ■ , Xk\, that is to say D is the smallest rectangle with sides 
parallel to coordinate axes, spanned by xi, . . . ,Xk- Furthermore, set Al to be the family 
of all vertical lines in [[D]] and for each I & Al — [[D]] let A*(/) be the topmost point of 
InD so that in particular A*(/i) = xi. With this notation, using Theorem 10 including the 
reversibihty statement there we see that the law of the conditional field A'^ is invariant and 
reversible for the following conditional version of the defective disagreement loop dynamics: 

• With intensity M.{dl)ds for I e A^ set ^s+ds '■— Is ® Should the proposed 
update fail or result in a configuration violating £[dl2.i X2] ■ ■ ■ ; dl^, Xfe] , keep ^s+ds — Is- 

• For each line I E Al extending an edge e of 7^ such that A*(/) G e, with intensity ds 
set js+ds '■= 7s0(^, A*(/)). Should the proposed update fail or result in a configuration 
violating S[dl2, X2; ■ ■ ■ •,dlk,Xk], keep js+ds = 7s- 

Note that the only difference compared to the DefDL dynamics in Subsection 6.4 is that 
we discard updates resulting in violation of S[dl2, X2; ■ ■ ■ ;dlk,Xk]. Observe also that 

• Under the so-defined directional decision rules the disagreement loops arising in the 
course of the dynamics are always single-branched (the second branch is chopped 
off by the boundary of D) and contain no cycles so that the updates in the above 
dynamics never fail (but note that they can nevertheless be discarded whenever 
violating £[dl2, X2] ■ ■ ■ ; Xjt]). 

• The only update of the dynamics which can lead to a transition from a configuration 
where £[dli,xi\ holds to a configuration where it does not hold is 0(/i,xi). Indeed, 
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in view of the directional decision rules, even if some disagreement loop annihilates 
a subsegment of li, this subsegment cannot reach Xi. 

• Likewise, the only update of the dynamics which can lead to a transition from a con- 
figuration where S[dli, Xi] does not hold to a configuration where it holds is (B{li, Xi). 

• For the update (B{h,xi) the only possibility to cause violation ofS[dl2, X2; ■ ■ ■ ;dlk,Xk], 
and thus to be discarded, is that the initial downwards creation phase growth of the 
disagreement loop from Xi continue with no directional updates along li until the 
loop reaches TZ[x2, ■ ■ ■ ,Xk]. Indeed, if a directional update occurcd prior to that, the 
disagreement loop would only unfold upwards and rightwards thereafter, never reach- 
ing TZ[x2, ... ,Xk] and consequently not violating S[dl2, X2; . . . ; dlk, Xk]- However, the 
probability that the initial segment of the disagreement loop for (B{xi,li) extends 
all the way down to 7l[x2, ... ,Xk] without directional updates is exp(— A^([[/i[xi 
7l[x2, ■ . . ,Xk]]]])) by the definition of the dynamics. Clearly, even if this happens, 
the violation of S[dl2, X2; ■ ■ ■ ; dlk, Xk] is not guaranteed, but possible. 

• Likewise, the update 0(/i,a;i), only possible on the event S[dli,Xi], may only cause 
violation of S[dl2,X2; . . . ; dlk,Xk] if the current configuration 7^ contains a segment 
along li joining xi to 7l[x2, . . . ,Xk]. On the event £[dli,xi] D £[dl2,X2; . . . ; dlk,Xk] 
this only happens with probability 0(exp(— 7W([[/i[a;i 7^[a;2, . . . ,2;^]]]]))) though, 
in analogy to the previous case. 

In view of these observations and by the definition of the dynamics we see that in its course: 

• On -iS[dli,Xi], transitions leading to S[dli,Xi] occur in s-time with intensity 

M{dk)ds (1 + (^e-M{[Mxi^n[x2,...,xm'j-^ . 

• On S[dh,xi], transitions leading to xi] occur in s-time with intensity 

ds{l + (e-M[[h[xi^nix2,...,xk]]]])-j-j _ 

Consequently, using the detailed balance we see that the probability of £[dli,Xi\ for A^, 
or equivalently the conditional probability of £[dli, Xi] given £[dl2, X2; ■ ■ ■ ;dlk,Xk\, is 

M{dh) (1 + (^Q-M([[h[x,^n^2,-,x,]]]])^^^^ 
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as required for (21). This completes the proof of the lemma. □ 

A repetitive use of Lemma 4 readily yields the following approximate factorisation 
theorem. 

Theorem 11 Assume the collection {h,Xi)i=i is such that Xi ^ TZ[xi^i, . . . ,Xk\ for all 
i and let Si :— Ai[[li[xi TZ[xi+i, . . . ,Xk]]]] if h hits TZ[xi+i, . . . ,Xk] and Si :— +oo 
otherwise. Then 



The purpose of the present section is to describe an alternative dynamics which can be used 
to simulate the empty-boundary field A^^f for /3 > 2. A particular feature of this dynamics 
is that it yields a perfect sampler, i.e. the sample it outputs comes exactly from the target 
distribution and not just its approximation. Moreover, for /3 large enough it extends to the 
whole plane thus allowing for a direct graphical construction (perfect simulation) of the 
whole plane thermodynamic limit A^'^. The idea underlying the graphical representation 
comes from Fernandez, Ferrari & Garcia [7, 8, 9] and in the context of homogeneous 
polygonal fields it has been developed in Schreiber [15, 16]. Here we present an extension 
of the construction to the general non-homogeneous setting. Our presentation splits into 
several steps. 

8.1 Free contour measure 

Choose an open bounded set D with piecewise smooth boundary and consider the space 
Cd consisting of all closed polygonal contours in D which do not touch the boundary dD. 
For a given finite collection (l) :— (Zi, In) of straight hues intersecting D denote by Cd{1) 
the family of those polygonal contours in Cd which belong to Td\${1) ■— ^d{1) H rD\$- We 
define the so-called free contour measure Qd on Cd by putting for C C. Cd measurable, 
say with respect to the Borel cr-field generated by the Hausdorff distance topology. 



a'^[dh,x,;...;dk,Xk] = llM{dLd J](l + 0(exp(-5,))). 




8 Contour birth and death dynamics 




e£CnCD(i) 



exp( 



M 




(22) 
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with Fin(L[D]) standing for the for the family of finite collections of lines intersecting D 
and where M* is the measure on Fin(L[£)]) given by dM.*{{li, ...,ln)) '■— nr=i-^(^^j)- 
Thus, in other words, 

QD{d9)^eM-L^id)) n M{dl[e]). 

eeEdges(6') 

For general /3 e R we consider the exponential modification of the free measure 
e^§{de) := exp(-[/? - l]L^{e))QD{de) = exp(-/5L^(^)) JJ M{dl[e]). 

eGEdges(9) (23) 

It is easily seen that the total mass 0^'(Cd) is always finite for bounded D for all G M. 
Indeed, taking into account that 1^(6) < card(Edges(6'))X([[i:>]]) we have in view of (23) 

oo 

k=o (24) 

Note that for all D the free contour measures ©d as defined in (22) arise as the respective 
restrictions to Cd of the same measure © on C := IJ^i ^(-n,n)2, in the sequel referred to as 
the infinite volume free contour measure. Indeed, this follows easily by the observation that 
Qdi restricted to coincides with for D2 C Di. In the same way we construct the 
infinite-volume exponentially modified measures ©[^], /? e M. It turns out that for /3 > 2 
a natural algorithmic random walk representation of ©[^^ can be given. This is done as 
follows. For each contour ^ e C let l[9] be its leftmost vertex, that is to say its extreme left 
point minimising the first coordinate, with possible ties broken in an arbitrary measurable 
way. Next, for x e denote by Cx the collection of contours 9 E C such that i[9] — x. For 
/3 > 2 let ©L'^' be the subprobability measure determined by the following construction of 
a Cj;- valued ©L^' -distributed random element 9 

• Simulate a continuous-time random walk {Zt)t>o gouverned by the following dynamics 

— choose two random lines Zi, Z2 meeting at x according to the distribution 

d[M X M]{{li,l2), /i n/2 e dx) 
2{{M)){dx) ' 

as in (7), 

— put Zq := and choose one of two possible initial directions along li with equal 
probabilities 1/2, 
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— between direction update events specified below move in a constant direction 
with speed 1, 

— while moving along the segment Z^t^t+dt] = ZtZt+dt-, update the movement direc- 
tion and start moving along a line / G [[Zft^t+^t]]] with probability 2M.[dl) and 
with equiprobable choices between the two possible directions along 

• Consider a killed modification of Zt by killing Zt 

— with intensity [(5 - 2]M{[[Z^t,t+dt]]]), 

— whenever Zf hits its past trajectory, 

• Draw an infinite loop- closing half-line /* along I2 beginning at x with equiprobable 
choices of two possible directions, 

• If the random walk hits the loop closing line 1* before getting killed, and the 
self-avoiding contour 9^. := ^^[Z'^"^]; /*] created by 1* and the trajectory of up 
to the moment of hitting /* satisfies i[6^] = x, then 

— with probability exp(— /3A^([[e*]])) output 9 :— 9^, where e* stands for the seg- 
ment of the loop closing line 1* from x to its intersection point with zf~'^\ 

— otherwise output 9 :— $. 

In all remaining cases set 9 

Note that outputting means failing in the above algorithm, and hence the resulting 
distribution is a sub-probability rather than probability measure. 

The following lemma provides a constructive representation for /5 > 2, and is close 
in spirit to Lemma 5.1 in [16] and Lemma 1 in [12]. 

Lemma 5 For P > 2 we have 

e[^] = 4 / e^f\{M)){dx). 
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Proof The directed nature of the random walk trajectories as constructed above re- 
quires considering for each contour 9 two oriented instances (clockwise) and 9*~ (anti- 
clockwise). By the construction of ©L^' the assertion of the lemma will follow as soon as 
we show that for each x e and 9 E Cx we have 

8((M))(dx)e-^^^(^*)p (zf-'l reaches I* and 9,[Z^'^-^hl*] e d9^) = e^'^\d9), 

(25) 

where e* stands for the last segment of 9'^ counting from x as the initial vertex, which is to 
coincide with the segment of the loop-closing line 1* joining to x its intersection point with 
zf~'^^; whereas 9^[Z^^~'^^;l*] is the self-avoiding contour created by I* and the trajectory 
of zf^'^'^ up to the moment of hitting as denoted in the construction of ©L^^ above. 
Indeed, the same relation holds then for 9'~, so adding versions of (25) for 9~' and 9^, 
which amounts to taking into account two possible directions in which the random walk 
can move along 9, will yield 2Q^^^(d9) on the right hand side, whence the assertion of the 
lemma will follow. 

To establish (25), we observe that the probability element 



P reaches I* and 9^[Z^-^^;l*] e d9^^ 



is exactly 



g^^^axp(-/,L-(»\e-))nA<(*N). (26) 

where ei, . . . , are all segments of 9 including e* — Ck, while l[ei\ stands for the straight 
line determined by e^. Indeed, 

• the prcfactor [8{{M)){dx)]-^ = [A[M x M]{{{h,l2) \ h n k E dx})]-^ times the 
product M.{dl[ei])A4{dl[ek]) comes from the choice of the lines /i and I2 containing 
respectively the initial segment of 9^ (counting from x) and the half-line as well 
as from the choice between two equiprobable directions on each of these lines, 

• for each of the remaining segments Cj, i = 2, . . . , A; — 1, the factor Ai{dl[ei]) comes 
from the directional update of the random walk Zt, with 2Ai{dl[ei]) due to the choice 
of the line and 1/2 due to the choice between two equiprobable directions along this 
line. 
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• finally, exp(— [/3 — 2]L-^(^ \ e*)) comes from kiUing in the course of the random walk 
whereas extra exp(— 2L-^(^ \ e*)) is due to the absence of directional updates along 
the segments of the walk, which yields exp(— /3L-^(^ \ e*)) when put together. 

Recalling (23) we see that the expression in (26) coincides with 

which yields the required relation (25). The proof is complete. □ 

To proceed, assume that the activity measure A4 admits a homogeneous upper and 
lower bound, that is to say 

C_/x <M< C+i^ (27) 

for some < C_ < C+ < oo and with /i standing for the standard isometry-invariant 
Haar-Lebesgue measure on the space [[K^]] of straight lines in R^. Recall that one possible 
construction of goes by identifying a straight line I with the pair (0, p) e [0, tt) x M, where 
(psin(0), pcos(0)) is the vector orthogonal to I, and joining it to the origin, and then by 
endowing the parameter space [0, tt) x R with the usual Lebesgue measure. Then, for each 
linear segment e we have by standard integral geometry 

2C_£(e) < M{[[e]]) < 2C+£(e) (28) 

because A*([[e]]) = 2i{e). Likewise, since ((//))(dx) = ndx as easily checked by straightfor- 
ward integration, 

C-ndx < {{M)){dx) < C+ndx. (29) 

Consequently, we conclude from Lemma 5 that under /3 > 2, the contour size exhibits 
exponentially decaying tails, which is a non-homogeneous counterpart of Lemma 1 in [15]. 

Lemma 6 Under the assumption (27), for j3 >2 we have 

Q^\{e] n Vertices(^) ^ 0, l{e) > R}) < 4C+7rexp(-2C_[/3 - 2]R)dx, 

where the event {lia; fl Vertices (^^) 7^ 0} is to be understood that a vertex of 9 falls into dx. 
Moreover, there exists a constant s > such that, for /3 > 2, 

e['^l({^; Oelnte, e{e) > R}) <e^p{-2C4f3 -2 + e]R + o{R)), 

with Int 9 standing for the region enclosed by the contour 9. 
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Proof The first assertion follows by the construction of the random walk Zt where, in 
view of (28), the kiUing intensity is at least 2C_ times the length element covered under 
the present assumptions, whereas the extreme left vertices of contours have their intensity 
bounded by ^{{M)){dx) < 4C+7r in view of the integral formula in Theorem 5 and of (29). 
To get the second assertion observe in addition that during each unit time of its evolution 
the random walk Zt has some uniformly non-zero chance of hitting its past trajectory, see 
also the proof of Lemma 1 in [15] . □ 

8.2 Polymer representation 

To proceed we let V^w] be the Poisson point process in Co with intensity measure O^^. 
It follows then directly by (23) in view of (2) and (3) that, for all (3 e R, coincides 
in distribution with the union of contours in Pc,[/3] conditioned on the event that they are 
disjoint so that 



where the conditioning makes sense because ©J5^ (Cd) is finite as shown in (24) above, see 
also [15], Section 2.2. In particular, in analogy to Subsection 2.2 and Theorem 2 ibidem, 
the law of A^^i^ is invariant and reversible with respect to the following contour birth and 
death dynamics (7^ ) onVo- 

(C : birth[/3]) With intensity Q^§{de)ds do 

• Choose a new contour 9, 

• If 6* n 7f = 0, accept 9 and set 7^^^^ := 7^ u 6*, 

• Otherwise reject 9 and keep 7^^^^ '■— 1^: 

(C : death[/5]) With intensity 1 • ds for each contour 9 e 7^^ remove 9 from 7^^ setting 



Moreover, jC{Aj} ) is the unique invariant distribution of the above dynamics, see Theorem 
2 in [15]. Clearly, both the representation (30) and the above dynamics can be regarded 
as a kind of polymer representation for A'^''^ . 




(30) 




'f \ 0. 
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Unlike the results for free contour measures, the polymer representation is valid for all 
values of (3. In this context it is helpful to note that, by (30), the probabihty P(^^0^ = 0) 
is not smaller than exp(— ©^'(Cd))- On the other hand, comparing with (2) shows that 
^(-^Dicf ^ ^) [^iy\^%]~'^ where, recall, Z^^^ is the partition function for ^;^0. Thus, by 

(24), 

Z^'^<exp(e['^](Cz.))<oo. (31) 

Likewise, 

Z^'^ < oo (32) 
which can be proven along the same lines, see also the proof of Corollary 2 in [15]. 

8.3 Graphical representation 

These observations place us within the framework of the general contour birth and death 
graphical construction as developed by Fernandez, Ferrari & Garcia [7, 8, 9] and as sketched 
below, see ibidem and [15] for further details. Choose (3 large enough, to be specified below. 
Define !F{C) to be the space of countable and locally finite collections of contours from C, 
with the local finiteness requirement meaning that at most a finite number of contours can 
hit a bounded subset of R^. On the time-space R x T{C) we construct the time-stationary 
free contour birth and death process {Qs)s&. with the birth intensity measure given by O'^l 
and with the death intensity 1. Note that free means here that every new-born contour 
is accepted regardless of whether it hits the union of already existing contours or not, 
moreover we admit negative time here, letting s range through M rather than just R+. 
Observe also that wc need the birth measure 0'*^] to be finite on the sets {dEC\6nA^ 
0} for all bounded Borel A C in order to have the process {qs)s&l well defined on 
M X J^{C). By Lemma 6 this is ensured whenever /? > 2. To proceed, for the free process 
{Qs)sm. we perform the following trimming procedure. Wc place a directed connection from 
each time-space instance of a contour showing up in {Qs)sm. and denoted by x [so,Si), 
with 9 standing for the contour and [so^si) for its lifespan, to all time-space contour 
instances 6' x [sq,s'^) with 6' \^ 6 ^ 0, Sq < sq and s'^ > sq. In other words, we connect 
6 X [so, Si) to its ancestors understood as those contour instances which may have affected 
the acceptance status of 6* x [sq, Si) in the constrained contour birth and death dynamics 
(C) as discussed above. These directed connections give rise to directed ancestor chains 



50 



of time-space contour instances, following Fernandez, Ferrari & Garcia [9] the union of 
all ancestor chains stemming from a given contour instance 9* — 9 x [so,Si), including 
the instance itself, is referred to as its clan of ancestors and is denoted by An(^*). More 
generally, for a bounded region U in the plane we write Ans(C/) for the union of ancestor 
clans of all contour instances 9 x [sq, Si) with 9 (lU $ and s e [sq, Si)- Lemma 6 allows 
us to apply the technique of domination by sub-critical branching processes, developed in 
[7, 8, 9], in order to conclude that there exists /3g such that for each (3 > (3g there exists 
c := c(/3) > such that 

P(diamAn,(B2(x, 1)) > R) < exp(-ci?), s e R, x e R^ (33) 

with M2{x, 1) standing for the ball of radius 1 in centred at x. In the sequel we shall 
always assume that (3 > (3g, that is to say that (3 is in the validity region of the graphical 
construction. We see that for P > Pg all the ancestor clans are a.s. finite and we can 
uniquely determine the acceptance status of all their members: contour instances with no 
ancestors are a.s. accepted, which automatically and uniquely determines the acceptance 
status of all the remaining members of the clan by recursive application of the inter-contour 
exclusion rule. In this case, discarding the unaccepted contour instances leaves us with a 
time-space representation of a time-stationary evolution (7s)sg]R on T{C\ which is easily 
checked to evolve according to the whole-plane version of the dynamics (C) above. In full 
analogy with Section 4 and Theorem 4 of [15] we see that for all s e IR the polygonal field 
75 coincides in distribution with the thermodynamic limit (see Section 3 ibidem) for A^''^ 
without infinite contours, which is unique (see Corollary 4 ibidem). For definiteness we 
put A^''^ :— 70- It should be also observed that for each s e R the free field Qs coincides in 
distribution with the Poisson contour process Vqw]. Since almost surely we have 7s C we 
get the stochastic domination of the contour ensemble A''^'^ by Vquj]. Moreover, using the 
exponential decay relation (33), in full analogy with Theorem 4 (iv,v) we readily establish 
the following /3-mixing statement for A^'^. 

CoroIIairy 5 For P > Pg there exists c > such that for each A,B C the law of 
the restriction A^^^ of the thermodynamic limit A'^'^ to A U B differs in variational 
distance from the product of the laws of the respective restrictions A'^'^, A^'^ by at most 
0(cover(74)cover(S) exp(— cdist(74, S))), where cover(X) stands for the minimum number 
of unit balls covering X C R^ whereas dist{A,B) :— mix^A,yeBdist{x,y). 
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Indeed, this is easily checked by noting that the considered variational distance is bounded 
above by the probabihty that either Ano(A) extends further than dist(A, B)/2 away from A 
or Ano(5) extends further than dist(A, B)/2 away from B which happens with probability 
at most 0(cover(74)cover(S) exp(— cdist(^, 5))) as required. Furthermore, in analogy to 
the discussion following Corollary 4 in [15] we can use the stochastic domination of ^^'^ 
by Vqip] and apply Lemma 6 to conclude that 

Corollciry 6 With probability 1 there is only finite contour nesting in A^'^i that is to say 
no infinite chains of nested contours are present. 

Indeed, this comes from the fact that the expected number of contours surrounding a given 
point is always finite by Lemma 6. In addition, for (3 large enough this expected number 
can be made arbitrarily small, which clearly guarantees the presence of long-range point-to- 
point correlations and spontaneous magnetisation with contours interpreted as separating 
phases of different signs, as follows by a standard Peierls-type argument. 

We also consider finite-volume versions of the above graphical construction, replacing 
the infinite- volume birth intensity measure ©^''^ with its finite- volume counterparts for 
bounded and open D with piecewise smooth boundary. Clearly, the graphical construction 
yields then a version of the finite- volume contour birth and death evolution (C). For each D 
denote by (7^ )sg]R the resulting finite- volume time-stationary process on the space J-{Cd) 
of finite contour configurations in D and write {Q^)se^ for the corresponding free process. 
It follows as in Theorem 2 in [15] that 7^ coincides in distribution with A^'^ for each 
s e M. Likewise, coincides in distribution with VaXn] . 

By representing the measures as the corresponding restrictions of ©[^1 we obtain a 
natural coupling of all the processes 7f , f?f , 7s and Qs on a common probability space. We 
shall also consider A^'^ coupled on the same probabihty space by putting A^'^ = 7,f . 
Under this coupling we have the following corollary stating exponential convergence of 
A-^'^ to A^'^ as D T R^. 

Corollciry 7 Assume that (5 is within the validity regime of the contour birth and death 
graphical construction. Then for x & D such that B{x,l) C. D we have 

P {[A^''%.i.,i) ^ [<1|B.(.,i)) < eM-c[dist{x,dD) - 1]), 
with c as in (33) and with \ standing for the restriction operator. 
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Indeed, this follows easily by (33) upon noting that the event {[^"^'^JiBaCa;,!) 7^ [«^d^''^]|B2(x,i)} 
only occurs when Ano(B2(a:, 1)) % D. 
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